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STUDY  OF  ANIMAL  SIGNALt  AND  NEURAL  PROCESSING, 
WITH  APPLICATIONS  TO  ADVANCED  SONAR  SYSTEMS 


1.  INTRODUCTION. 


This  progress  report  summarizes  the  work  done  under  contract  No. 

N001 23-73- C- 1144,  The  report  extends  the  results  Obtained  under  two  previous  con- 
2>  3 

tracts.  AH  three  contracts  were  performed  under  the  auspices  of  the  Biosystem 

* 

Research  Department  of  the  Naval  Undersea  Center,  San  Diego,  California. 

The  report  begins  with  a  summary  of  past  results,  along  with  some  re¬ 
marks  that  pertain  to  newer  discoveries.  An  important  new  approach  to  the  analysis 
of  animal  echolocation  behavior  is  presented  in  Section  3.  This  approach  is  exciting 
because  it  relates  the  sonar  signals  used  by  bats  and  cetaceans  to  well-accepted  models 
of  the  mammalian  auditory  processor.  It  is  analogous  to  approaching  the  problem  from 
"away  around  on  the  far  side",'  a  vantage  point  from  which  animal  echolocation  behavior 
and  the  theory  that  has  been  used  to  explain  it  became  much  more  obvious  and  easily 
accepted. 


The  fourth  and  fifth  sections  consider  likelihood  ratio  testing  and  signal 
to  Interference  maximization  as  applied  to  our  model  of  the  animal  sonar  system.  It 
is  demonstrated  that  wideband  signals  constitute  a  more  practical  set  of  basis  func¬ 
tions  than  sinusoids,  when  the  receiver  is  completely  linear. 


*The  advice,  encouragement,  and  criticism  of  W.  E.  Evans  and  C,  S,  Johnson  (both  of 
NUC)  has  been  especially  helpful. 


1.  — Continued. 

The  sixth  section  describes  some  new  results  concerning  bet  signals.  It 
3 

was  previously  thought  that  bat  signals  are  used  to  convey  amplitude  modulation  {AM) 
information.  More  recent  results  indicate  that  bat  signals,  like  those  of  cetaceans, 
can  be  explained  in  terms  of  linear  filter  characterization.  This  finding,  however, 
does  not  negate  the  AM  hypothesis.  For  certain  carrier  signals,  it  happens  flat  AM 
can  be  described  as  a  linear  filtering  operation. 

The  seventh  section  describes  an  experiment  that  is  presently  being  under¬ 
taken  as  a  practical  test  of  the  theory. 

The  eighth  section  describes  two  alternative  models  of  the  auditory  system. 
Both  models  involve  a  set  of  constant-Q  filters,  but  the  model  suggested  by  the  author 
incorporates  filters  with  narrower  bandwidlhs,  which  are  viewed  as  the  basic  components 
from  which  the  constant-Q  filters  are  constructed. 

The  ninth  section  deals  with  the  problem  of  minimizing  the  volume  under 
the  wideband  ambiguity  function.  Volume  minimization  is  of  only  incidental  importance 
for  animal  echolocation,  but  is  of  interest  to  designers  of  conventional  sonars  when 
planar  targets  are  assumed. 

It  has  been  mentioned  that,  for  certain  carrier  signals,  AM  can  be  described 
as  a  linear  filtering  operation.  This  set  of  carrier  signals  is  composed  of  all  waveforms 
that  can  be  expressed  ftf?  finite  order  polynomials  over  a  bounded  frequency  interval 
[  0,  W  3,  The  advantages  of  these  waveforms  for  characterizatton  of  time  varying 
linear  systems  and  for  communication  are  discussed  in  Section  10.  Taylor  series 


-4 

i* 


--Continued* 


characterization  of  signal  spectra  provides,  a  viewpoint  from  which  one  can  easily  dis¬ 
cern  the  quantitative  equivalence  between  time  varying  filters  (or  targets,  or  channels) 
and  time  invariant  ones. 

The  mathematical  processes  involved  in  animal  echcloeation  are  intimately 
related  to  the  processing  of  speech  by  the  human  auditory'  system.  This  relation  is 
pursuedinSectlon  11,  where  a  new  method  of  extracting  information  from  speech  is 
suggested.  The  method  can  be  interpreted  as  a  speech  decoding  algorithm.  Such  an 
algorithm  can  be  applied  to  data  reduction  for  transmission  of  speech,  to  computer 
understanding  of  speech,  and  to  phoneme  synthesis. 

The  Appendix  includes  some  mathematical  properties  of  the  signals  that 
have  been  derived  to  explain  animal  eeholocation  behavior. 


Nations  d*  Sire  :  Motivations. 


Although  the  preceding  paragraphs  preview  the  contents  of  lids  report, 
they  may  not  provide  sufficient  reason  to  proceed  unless  a  particular  subject  has  caught 
the  reader’s  eye,  I*  is  advisable,  then,  to  consider  the  report,  from  the  viewpoint  of 
motivation. 

Ttie  U,  S.  Navy  has  spent  more  than  a  few  dollsro  tc*  support  this  research, 
and  the  author  has  spent  more  than  a  few  hours  to  carry  it  out.  Tc  what  ends  have  this 
time  and  money  bsen  spent?  The  report  gives  some  results,  but  what  were  the  questions 
and  concerns  that  motivated  tee  research?  Some  oi  the  question?  are  well  defined;  otters 
are  somewhat  nebulous.  The  questions  seem  to  fit  three  major  categories: 

L  Sonar  Systems. 

A.  Specific  Questions. 

How  should  a  sonar  system  deal  with  the  commonly  observed 
situation  in  which  a  target  echo  only  faintly  resembles  a  wideband  trans¬ 
mitted  signal?  Is  there  a  general  mathematical  description  for  such  a 
phenomenon?  Can  the  phenomenon  be  exploited  for  target  recognition 
or  clutter  suppression?  Can  it  be  circumvented  so  that  matched  filter¬ 
ing  can  be  applied  to  wideband  sonar  echoes?  What  signals  and  filters 
are  optimum  for  (a)  exploiting  or  (b)  circumventing  the  phenomenon? 

B,  Nebulous  Questions. 

Can  an  automatic  man  mate  sonar  system  duplicate  {or  better) 
toe  target  recognition  and  clutter  suppress!®  capability  of  cetaceans 
and  bats?  How  would  such  a  system  work? 


i- 


Animal  Echolocatieo, 


specific  Questions, 

is  there  a  restricted  set  of  signals  that  encompasses  most 
animal  echolocatton  waveforms?  Is  there  a  general  theory  that  is 
applicable  to  both  porpoise  and  bat  sonar?  What  is  the  meaning  of 
existing  neurophysiological  and  psychoacoustic  data,  in  terms  of  a 
mathematical  model  of  the  animals*  receivers? 


Nebulous  Questions, 

If  animal  eeholoeation  behavior  is  mathematically  explained, 
what  then?  Will  neurophysiologists  and  paychoaeousticians  be  given 
a  new  set  of  interesting  experiments  to  perform,  in  order  to  discover 
how  the  theory  is  implemented?  Can  neurophysiological  data  reveal 
new  processing  implementations  that  would  be  useful  to  the  engineer? 
Does  an  uniterstanding  of  animal  echoiocation  lear  to  an  understand¬ 
ing  of  other  auditory  processes,  such  as  speech  recognition? 


Communication  Theory,  Medicine,  and  Otter  Applications, 


Specific  Questions, 


How  do  the  answers  to  the  above  questions  directly  benefit  the 
analysis  mid  synthesis  of  man  made  communication  systems,  as  well 
as  sonar  systems?  Are  there  some  straightforward  applications  to 
medicine,  e*g,,  to  diagnostic  ultrasound? 
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1*1  — Continued. 

B.  Nebulous  Questions. 

How  can  other  branches  of  engineering  or  science  benefit 
from  an  understanding  of  animal  echolocation  systems? 


Most  of  the  above  questions  are  answered*  with  varying  degrees  of  con¬ 
fidence  and  completeness,  in  this  report*  It  is  the  author's  hope  that  the  report  art 
only  answers  some  past  questions,  hit  generates  insights  for  interesting  and  useful 
Mure  work.  It  is  also  hoped  that  the  reader  now  has  sufficient  motivation  to  investi¬ 
gate  this  report*  as  well  as  related  research,  in  further  detail. 


"The  prospects  of  using  animal  echolocation  are  inexhaustible  in  the  area 

of  technology,  engineering  development  and  in  all  those  problems  which  have  been 

raised  by  bionics.  Numerous  articles  and  bodes  have  teen  written  on  this*  and  diverse 

judgments  have  been  offered  at  special  conferences  and  symposia.  But  up  lo  now  they 

have  reflected  usually  hypotheses,  searches,  desires,  and  dreams  rather  than  specific 

results*  We  might  recall  here  the  wise  words  of  the  professor  at  Leningrad  University, 

A,  A.  Ukhtomskiy:  ’Man  is  a  very  powerful  being:  if  he  begins  to  dream  seriously, 

31 

tBen  this  means  that  sooner  or  later  the  dream  will  come  true'. " 


ESL-PR115 


2.  REVIEW  OF  PREVIOUS  RESULTS. 


2.  1  A  Signal  Processing  ModeL 


Tbe  cetacean  sonar  signal  processor  Is  modelled  as  shown  in  Figure  2-1. 
Auditor}'  threshold  data  (Figure  2-2}  seems  to  imply  the  existence  of  the  whitening  filter 
shown  In  the  figure  2-1.  For  a  constant  internal  threshold  T  ,  the  externally  measured 
threshold  of  tearing  Tp  (cu)  should  satisfy  the  equation 

T£  <«)  W(w)  Tj  (2-1) 

where  W  (u>)  Is  the  transfer  function  of  the  whitening  filter.  Given  the  shape  of  W  (&?) 

from  Figure  2-2  and  Equation  (2-1),  one  would  expect  W  (a?)  to  vary  Inversely  with 

S_,  (u?),  tte  ambient  noise  power  spectral  density.  Figure  2-3  indicates  that  W  (&>) 

N 

indeed  behaves  as  a  whitening  filter* 


2. 2  _ Description  of  Objects;  Target  Transfer  Function. 


Many  engineers  have  suggested  that  radar/  sonar  echoes  can  be  described 

4_g 

as  linear  transformations  of  transmitted  waveforms.  ‘  This  description  Implies  that 
tte  Fourier  transform  of  tte  echo,  E  (w),  can  be  expressed  as  the  preduct  of  the  target 
transfer  function,  T  (w),  and  tte  Fourier  transform  of  tte  transmitted  signal,  U  <w); 

EteO  =  T  (cj)  U(w>  .  (2-m 

Simple  geometrical  forms  (spheres,  cylinders,  cones,  etc. )  have  radar 
2 

or  sonar  cross  sections,  J  T  (w)  j  ,  that  can  generally  be  represented  by  only  a  few 


fi*«»  2-1,  The  Optimum  Linear  Receiver  lor  uncorrolated  signals,  u 
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2.  2  —  Continued, 


The  first  set  consists  of  all  signals  that  satisfy  condition  (2-5).  The  second  3et  consists 

2  3 

of  signals  that  ate  sensitive  to  scale  mismatches  of  order  k,  k  ,  k  ,  etc, 

3 

It  can  be  shown  that  the  whitening  filter  in  Figure  2-1  has  no  effect  upon 
target  characterization  with  the  signals  (2-8),  provided 

W  (w)  «  w 2  (2-9) 

as  shown  in  Figure  2-3.  Each  filter  in  Figure  1  is  then  matched  to  a  differently  scaled 
version  of  the  transmitted  signal,  and  the  output  of  each  filter  estimates  the  value  of  a 
different  target  coefficient  in  Equation  (2-3) . 


Tursiopa  Echolocation  Signals 


The  measured  waveforms  *  in  Figure  2-6  were  recorded  while  a  Tursiops, 

vision  occluded  by  opaque  rubber  cups,  was  discriminating  between  two  acoustically 
12 

different  targets. 

The  theoretical  signals  in  Figure  2-6  correspond  to  different  versions  of 
(2-8),  as  listed  in  Table  2-1.  The  animal  whose  signals  are  shown  in  [Figure  2-6  was 
torn  wild,  and  demonstrated  a  well  developed  ability  to  discriminate  between  targets 


*  All  porpoise  signals  shown  in  this  report  were  measured  by  W.  E,  Evans  of  NUC, 
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Table  l"  1 .  Maths  mat  tea’  Expressions  for  the  Theoretical  Signals  in  Figure  2-6. 

1.  U(oj)  =  a°gU)Z  [l  +  f.lojje'38^10^] 

.  Cl-(.0Sj)e‘J4"logw]  eJ-(-2»Mg«+»/8) 

2.  U  (w)  =  oT°* 9  e“(log  w>2  [I-  (0. 2j)  e~3  (8“  log  w  '  7/43  ] 

-  i  2#  log  u 

3.  U(U)  =  u"11-2  e'4  flog  “)2  e-HHogu 

-0.  2  sin  (8  S’  log  «;) 

4.  U  («)  =  aT5*4  g-2^0^^)  e-j6ff  logo> 

-0,1  cos  (4ir  log  U’>  -j  rr/2 

•  8 

5.  U(w)  =  w“11,2  e~4  0°S  w)2  e-j4-  Iog^ 

-0.2  sin  (Sir  log  t)  -  j  jr/2 

•  v 
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% 

by  means  of  echolocation.  At  the  time  the  signals  were  recorded  this  ability  had  become 
perfected  even  further*  so  that  comparatively  few  pulses  were  used  to  perform  the  dis¬ 
crimination. 

The  measured  waveforms  in  Figure  2-7  were  recorded  when  a  dolphin  was 
performing  range  discrimination  between  two  identical  targets.  The  animal  in  this  case 
was  a  Tursiops  that  had  been  born  and  raised  in  captivity.  When  this  inexperienced 
animal  was  first  induced  to  accept  the  rubber  cups  that  occlude  vision,  he  avoided  large 
objects  but  was  very  cautious  about  getting  too  close  to  them.  The  animal  could  not 
initially  find  or  recognize  dead  fish  and  exhibited  an  obviously  startled  reaction  when 
he  happened  to  bump  into  one.  The  dolphin  learned  quickly,  however,  and  could  soon 
recognize  a  fish  or  a  hand  in  the  water,  U  was  interesting  that  a  fish  held  in  a  hand 
constituted  a  new  recognition  problem,  even  though  fish  and  hand  could  be  recognized 
separately. 

Even  after  developing  the  ability  to  discriminate  between  targets,  the  inex¬ 
perienced  animal’s  sipials  were  occasionally  unorthodox,  as  shown  in  Figure  2-7.  Although 

the  two  measured  signals  in  Figure  2-7  do  rot  appear  similar  to  those  in  Figure  2-6, 
the  theoretical  waveforms  can  still  be  expressed  in  terms  of  its  functions  {2-8).  The 
parameters  of  the  theoretical  signals  are  given  in  Table  2-2. 

The  theoretical  waveforms  in  Figures  2-6  and  2-7  are  very  similar  to 
those  used  by  Tursiops.  This  similarity  suggests  that  the  porpoise  signals  are  "good 
waveforms  for  target  description, "  as  depicted  in  Figure  2-5.  The  theoretical  wave¬ 
forms  that  coincide  with  Tursiops  sisals  can  be  further  analyzed  to  determine  their 


e~j  [  S-  {.  94>  log  sa  +  .4 *  1 

.7  q-.  94  (log  co)2  [l*2.0e'3  L4r(-^ 

ri-.32Se-i[M(-94),°S“”/43 


Figuro  2-7,  Theoretical  and  Measured  Tara  lops  Sonar  Signals 

(Inoxparloncad  Animal  Born  in  Captivity) 
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efficacy  for  target  description  and  clutter  suppression,  as  well  as  their  resistance 
to  confusion  when  objects  are  moving. 

The  wideband  ambiguity  function  of  the  first  theoretical  signal  in  Figure 
2-0  is  shown  in  Figure  2-S.  Figure  2-8  shows  the  region  of  the  stretch  axis  where 
typical  target  velocities  {e20  knots)  would  appear.  The  response  of  the  filter  U*  (ai) 
to  energy  normalized  echoes  of  the  form  k~*  U  (u/k  )  for  n  =  1  and  n  =  2  cor¬ 
respond  to  profiles  drawn  along  the  dotted  lines  in  the  figure. 

The  ambiguity  function  can  be  interpreted  as  the  response  of  the  filter 

U*  (&) ,  matched  to  the  transmitted  signal,  to  energy  normalised  echoes  that  are 

scaled  versions  of  the  transmitted  signal.  With  this  interpretation  it  is  evident  that  the 

signal  will  work  well  for  detecting  large  planar  objects  surrounded  by  small  Rayleigh 

2  18 

seatterers.  The  sonar  cross  section,  j  T  (a?)  |  *  of  a  Rayleigh  scatters r  varies 

4 

as  ta  ,  Since  such  a  scatterer  is  small  in  comparison  with  a  wavelength,  T  («)  will 

have  negligible  phase  variation  as  a  function  of  frequency.  Therefore  T  to)  Is  pro- 
2 

portional  to  ~  .  Tie  fiiieris  response  to  a  planar  target  Is  at  s  =  1  in  Figure  2-8, 

o 

while  the  response  to  a  Rayleigh  scatterer  is  at  s  =  k  ", 

The  ambiguity  function  can  also  be  interpreted  as  the  response  of  the 

^  n /  2  « 

filter  k  U*  tw/k  i  to  other  seated,  energy  versions  of  D  to).  If  cat* 

is  looking  for  a  small  scatterer  against  a  large,  flat  reflecting  surface,  the  filter 

*  n 

function  is  k  “  U*  (£c/k*)  and  the  clutter  is  the  unsealed  signal  with  transform  U  (w). 

The  response  of  this  filter  to  an  energy  normalized  echo  from  the  flat  reflector  is 


Poult  I  on  n  of  k1  and  kz  Profllou 


lasn  gwen  by  the  s  =  ic  profile  to  Figure  3-g,  wmle  the  itiparise  to  the  small  t 
terer  Is  given  by  tbe  s  =  1  profile. 


The  ambiguity  function  of  the  theoretical 


the  waveform  should  work  well  for  detecting  a  small  object  amid  planar  clutter.  1 
Tursiops  appear  capable  of  locating  and  identifying  a  small  target  that  Is  position* 
against  a  highly  reflective  surface,  e.g..  fish  floating  at  air/ water  Interface,  fis] 
in  contact  with  a  rock,  or  fish  on  a  fiat  bottom. 


The  shaded  region  of  the  atadboity  plane  in  Figure  3-8  illustrates  tha 
little  confusion  is  introduced  by  moving  facets,  since  tbs  sews  factor  associated 


Is  ssmefently  scale  sensitive  for 


but  lacks  the  timi 


taudwldtfc 


for  resolv 
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2. 3  —Continued. 


functions  are  reasonably  well  icnown  (or  If  they  can  be  measured)  then  it  should  be  pos¬ 
sible  to  discriminate  against  interfering  reflectors  even  when  the  reflectors  have  range 
and  velocity  identical  to  those  of  the  target. 

Let  V*  (co)  be  the  transfer  function  of  the  filter  that  is  used  to  process 
received  echoes,  and  let 


Ec  (w)  =  C  (co)  U  (co) 


(2-11) 


be  the  echo  from  a  typical  clutter  reflector.  If  v  (t)  is  the  inverse  Fourier  transform 
of  V  (co),  then  the  response  of  the  filter  to  a  given  clutter  echo  is 


R(r) 


v  (t)  ec*(t  +  r)  dt 


(2-12) 


where  e  (t)  is  the  inverse  Fouri  Jr  transform  of  E  (co). 
c  c 


In  terms  of  spectra. 


R(t) 


(1/  2j r) 


f  V  (co)  Ec*  (co)  e~iWT  dco 


V  (co)  IF  (co)  C*  (co)  e'JW  r  dco 


(2-13) 
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Let  p  (r)  be  a  probability  density  function  that  describes  the  range  distribution  of  clutter 
reflectors,  where  r  is  measured  relative  to  the  target  The  expected  output  power 
(clutter  response)  of  the  filter  is  then 


clutter  response  =  E 


P  (t)  1  R  <r)  |  dr 


(2-14) 


where  E  is  the  total  energy  of  clutter  returns.  Since  r  is  measured  relative  to  the 
c 

target,  r  =  0  corresponds  to  the  arrival  time  of  the  target  echo.  The  response  of  the 
filter  to  the  target  echo  is 


target  response  =  E  j  (1/  2?r) 
9 


where  is  the  signal  energy. 


/» 

T  (£0)  U  (w)  V*  (w)  J  w  T  dw  I  2, 


(2-15) 


Since  the  maximum  response  is  assumed  to  occur  at  r  =  0, 


maximum  target  response  = 


Es  |  <1/2*} 


T(&?)U(w)V*(u?)d&i  |  “ . 


(2-16; 


I m 


Assuming  that  the  noise  is  approximately  white  over  the  receiver’s  bandwidth,  the 
expected  noise  response  of  the  filter  is 
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— Continued, 


noise  response  =  <Nq/  2)  '  (1/  2  «; )  j  V 


(w)  j  2  dee 


(2-1?) 


where  N  /  2  is  the  noise  power  spectral  density, 
o 


Interference  is  the  sura  of  the  noise  response  and  the  clutter  response. 
The  signal-to- interference  ratio  is  then 


Es  laT 


r 

J 


T  (o>)  U  (w)  V*(«)  do? 


SIR  = 


N  1  f  jV(w>j2dw  +  Ec  l  p(r)|  R(r)  j  “dr 

n  *  n*.  ** _ 


2  fr  -® 


JiL 

N 


i: 


T(a)  U  (cc)  V*  (it})  dec 


/  11 


2E  ” 

V  (y)  |  dec  +  ~ j  I  V  (x)  H  (y,  x)V*(y)  dxdcc 
°  (2- IS) 


where 


M  = 


^  r 


S 

I 


T  (y>  U  (y)  v*  (w)  dy 


2-24 
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-Continued. 


and 


H  fatx) 


f 


p  (r)  V  fa,  C  fa  V*  (x)  C*  (x)  e  -  <x"^Td r 


(2-20) 


The  filter  function  V  (a?)  that  maximizes  SIR  will  maximize  the  numeratci* 

of  (2-1S)  while  constraining  the  denominator  to  he  small.  For  energy  normalised  filter 

21 

functions,  this  approach  leads  to  a  constrained  version  of  the  Schwarz  inequality. 

The  solution  is 

fa 

V  fa  =  <2p/N  )  T  fa  U  fa)  -  (2E  /  2  ff  N  )  j  V  <x)  H  (&>,x)  dx. 

O  C  O  I 

(2-21) 


If  clutter  is  assumed  to  be  uniformly  distributed  in  range,  then 


U  fa,  x)  ~  2~  p  (0)  U  fa)  C  fa)  U*  (x)  C*  (x)  6  (x  -  (a) 


(2-22) 


and 


V*  fa) 


ii  T*  fa)  V*  fa) _ 

E  p  (0)  j  U  fa)  C  fa)  j  “ 


(2-23) 


Using  (2-10)  and  (2-11).  the  maximum  SIR  receiver  for  white  noise  can 


be  written 
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Et*  M 

F(o>)  = - - $ 

kr  (Ec  (a*)  I 


(2-24) 


Equations  (2-23)  and  (2-24)  give  the  best  filter  function  for  a  given  signal. 
The  next  step  is  to  find  the  best  signal  for  a  given  filter. 


Since  we  are  working  with  unit  energy  functions. 


I 


{  V  (o>)  }  dw  = 

Z.&  "1-S? 


CD 

J  IUWI2  a»  -  „  2r.  (3-25) 


Equation  (2-18)  can  than  be  rewritten 


"  / 

O  A, 


T  ha)  V*  ((a)  U  (o>)  dvJ 


SIR  = 


I  U(cu)|2dW+  iTf~  | 

°  i 


i  U*  (x)  G  {(a.,  x)  U  (a,)  dxdtu 

(2-26) 


P  IT)  V*(w)  C  (a)  V  (x)  C*  {%)  e 


■j  (x-tu)  r 


G  (o,x)  s 


(2-?/) 


--Continued, 


£r*  %t 


Comparison  of  Equations  (2-26,  2-2?)  with  Equations  (2-18,  2-20)  reveals  that  when 

V  (w)  is  replaced  by  U*  («)  (2-28) 

and 


U  (tu)  is  replaced  byV*  (no)  (2-2S) 

the  SIR  is  unchanged.  The  substitutions  (2-28,  2-29)  can  then  be  applied  to  Equation 
(2-23),  yielding 


U  (cu) 


_ n  T*  (**)  V  (o>) _ 

!*  2 
~ -4  Ep(0)  |  V  (u)  C  (W)  | 


(2-30) 


Indeed,  the  expression  (2-26)  is  maximized  when 


U*  (a>) 


(2p/No)  T  (w)  V*  (u)  -  (2Ee/2  «  KM 


U*  (X)  G  (u,  x)  dx, 

(2-31) 


and  whan  p  (r)  is  uniformly  distributed  in  range,  Equation  (2-31)  becomes  Identical  with 
Equation  (2-30), 
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Given  an  initial  signal  with  Fourier  transform  Ut  {(a),  the  best  filter  func¬ 
tion  V  *  (a?)  is  given  by  Equation  (2-23),  The  best  signal  U  (cu)  corresponding  ,*o  V  (tu) 
1  «  1 

is  then  given  by  Equation  (2-30),  One  can  thus  find  the  optimum  signal-filter  pair  by 
iterative  substitution  into  Equations  (2-23)  and  (2-30),  The  optimum  combination  will 
result  in  negligible  change  when  additional  iterations  are  applied  in  an  attempt  to  obtain 
better  signal  and  filter  functions. 


The  above  procedure  gives  an  optimum  signal-filter  pair  for  clutter  sup¬ 
pression  on  the  basis  of  reflector  shapes,  for  a  white  noise  background.  The  usual 
background  encountered  in  the  open  ocean  and  in  harbors  is  generally  colored,  however, 
as  shown  In  Figure  2-3.  For  a  noise  spectrum  S„  (cu)  that  has  a  colored  comnonent 

S  (u).  the  optimum  filter  function  changes.  Let 
c 


SM<W>  V2  +  s  * 

n  o  e 


(2-32) 


so  that  the  noise  response  in  (2-17)  becomes 


noise  response  =  (N^/2)  *  (1/ 2~)  I  |  V  (cc)  | “  dec 


i 


+  (1/2-) 


S  M  |  V(w)  I  “  dec 


(2-33) 


2-28 


9 


— Continued, 


When  the  extra  term  In  (2-33)  Is  Included  in  the  denominator  of  the  SR  in  (2-18), 

the  filter  function  that  maximizes  SIR  must  satisfy  the  equation 


v  m 


r  (a)  u  (<u)  - 


2E 

_ _ c 

2s-N 


o 


v  (x)  H  (03,  x)  dx 


-  (2/Ny  S c(03)  V  (03)  , 

i.e,,  an  extra  term  is  aided  to  the  right  land  side  of  Equation  (2-21), 


(2-34) 


Again  aiming  that  clutter  Is  uniformly  distributed  In  range,  the  optimum 
filter  function  becomes 


_ _ I*  T*  fa)  U*  fa) 

N ‘  ”  - - - 

~T  +  sc  <a,)  -  Ecp  {0)  1  u  <w>  c  (»)  i 


(2-35) 


Since  the  colored  mdse  terra  In  (2-33)  is  independent  of  the  signal,  the 
expression  (2-30)  for  the  optimum  signal  is  unchanged.  Nevertheless,  a  simultaneous 
solution  of  Equalises  (2-30)  aid  (2-35)  for  an  optimum  signal-filter  pair  will  generally 
result  In  a  different  signal  function  wteu  colored  noi so  is  introduced. 

As  In  the  clutter-free  situation,  it  appears  that  the  optimum  filter  Is 
matched  to  the  signal  when  noise  is  wWte,  but  a  mismatched  filter  is  optimum  when 


noise  is  em  «Mtf,  Tola  obsefvsfiesi  Is 


Uf 


MeMteaL 


asumoMss,  fcoi;;  i^WCi  sn-d 


ir  that  V  ft 


ttsoss 


lejl  are 


9  +  6CP@)  !  v  c»j  c  ft 


P  T  *  ft?) 


Taking  the  magnitude  of  both  sides. 


j  V  ft?)  r  =  |  U  (a?) 


i!  T(«)  !-No/2 
ecp(o>  jcft?)  r 


12-371 


In  frequency  intervals  where  N^/  2  >p  |  T  ft?)  { ,  ns  useful  information  ean  be  obtained 
for  making  a  decision  about  the  presence  or  absence  of  a  target  Therefore, 


2  ,  max  f  f»  |  T  (w)  j  -  N./  2 ,  0  j 

i  v  ft?)  i  =  j  u  {«)  r  =  — - — 

e_  p  m  i  c  ft)  r 


i'2-3S) 


If  V  fa?)  and  U  ft?)  are  assumed  to  be  identical  for  toe  colored  noise  situa¬ 
tion,  the  above  procedure  leads  to  a  contradiction;  (2-30)  and(2-35)  give 


a  j  T  (or)  J  =  ~  +  S  ft?)  +  E_  p  (0)  [  V  ft?)  C  ft) 


-f  +  E  p  (0)  i  v  ft?)  C  ft?)  i 

2C  • 


*2-39) 


two  quantities  on  the  right  hand  side  of  (2-39)  cannot  be  equal  unless  S  to)  -  0. 
ttched  filter  is  only  optimum  wten  the  additive  noise  has  constant  power  soeetral 
ity*  For  colored  noise,  a  mismatched  filter  results  In  a  larger  signal  to  Inter- 


1 

I 

I 
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3. _ LINEAR  FILTERS  FOR  TAYLOR  SERIES  SIGNAL  ANALYSIS,  WITH 

APPLICATIONS  TO  THEORIES  OF  HEARING  AND  ANIMAL  ECHO- 
IOCATTON. 


How  does  one  determine  the  Taylor  series  spectral  coefficients  of  any 
given  time  function  ?  The  method  that  first  springs  to  mind  is  to  evaluate  the  deriva¬ 
tives  of  the  function's  Fourier  transform  at  to  =  0,  This  method,  however,  las  certain 
disadvantages  which  will  be  discussed. 

A  better  scheme  would  be  to  find  a  linear  filtering  operation  such  that  each 
coefficient  Is  determined  by  measuring  a  maximum  filter  response.  Such  an  operation 
can  indeed  be  constructed.  The  restfiling  processor  is  a  sank  of  constant  Q  filters. 


The  Importance  of  this  filtering  technique  is  that  the  resulting  processor 
corresponds  closely  to  psvchoaceustic  models  of  the  raammalian  auditory  spectrum 
analyzer.  The  correspondence  implies  that  the  sound  processors  of  porpoises,  bats, 
and  people  are  all  designed  for  Taylor  series  spectral  analysis.  The  efficient  utili¬ 
zation  of  these  processors  as  sonar  receivers  involves  the  transmission  of  a  particular 
set  of  echoic-cation  si  goals.  This  signal  set  is  the  same  as  the  me  described  in  Sec¬ 
tion  2. 


la  summary,  one  can  start  from  the  premiss  that  signals  are  to  be  analysed 
ta  terms  of  their  Taylor  series  spectral  coefficients.  With  this  premise,  one  can 
derive  both  the  signals  and  the  sonar  processors  that  are  employed  for  animal  echo- 


location. 
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3. 1  Taylor  Series  Signal  Analysis. 

Sup  pc  sc*  ft  receded  signal  f  ft)  is  to  be  described  in  terms  of  the  complex 

poser  series  coeffictats  of  its  Fourier  transform  F  fa?) ,  i.  e.  f  in  terms  off  ,  f  , 

o  1 

. . .  }  f  f  . . .  ^  -a  here 

3 


F  M  = 


53 

2 

3=0 


(£*D 


A  s:  raightfoywa  rd  way  to  obtain  tea  Taylor  coefficients  ntt  Is  to 

B  '  S 

multiply  f  ft)  by  (j  t)  and  measure  the  d.  c.  level  of  the  resulting  function,  S.  e.  *  to 

evaluate 


d 


m 


F  (w) 


m  I  f  .  #3-2) 

m  ‘  - 


Such  a  procedure  is  dlsads'antageo-as  because  It  involves  a  ncclioear  oper¬ 
ation  is  the  time  domain  mid  because  It  assumes  that  the  time  of  arrival  of  f  (t)  Is 
known.  B  is  also  only  reasonable  if  f  (t)  tats  a  known  duratioa  %  so  that  {3-2}  can 
bo  determined  by  the-  operation 

T 

(Jt)*“fp)dt  =  m  J  f  r  ,  (3-3) 


1-2 


The  interval  [0,  T  1  over  whfeh  f  (t)  is  not  identically  zero  most  be  weQ 
denned  a  priori,  unless  there  is  no  extraneous  wise. 

When  white  noise  is  added,  the  ratio  of  signal  power  to  noise  power  at  the 
output  of  the  processor  (3-3)  Is 


o 


where  N_/ 2  is  the  power  spectral  density  of  the  added  noise.  For  T  =  1,  (3-4)  indicates 
that  the  ignal  to  noise  ratio  Increases  with  na,  so  that  the  higher  order  components  are 
more  accurately  estimated  than  the  lower  order  ones.  Tie  lower  orcer  terms,  however 
are  generally  more  Important  than  the  higher  coefficients  if  F  (as)  Is  to  be  approximated 
over  a  bounded  frequency  Interval  with  a  finite  number  of  coefficients.  The  dependence 
of  signal  to  noise  ratio  o*t  the  order  of  tbs  estimated  coefficients  must  be  listed  as  a 
disadvantage  of  the  method  described  by  (3-3). 

The  above  disadvantages  will  be  overcome  If  we  caa  Had  a  set  of  uncor- 
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To  say  that  the  unit  ener gy  functions  d  (_■'  are  mtorrelsted  Is  to 

B 


If  *  isi  2 

“7*  I  #{&?)♦  (wi  da*  |  « 

2  -  J  n  m  * 


P-Q 


Condition  (3-6)  is  somewhat  stronger  than  the  us *12!  qassl-orthogfiiialtty  csr.dltu 


.1 


d  tel  d  te)dw  =  {  a«mter«  l» 

"  “  mfu. 


The  iwerrsls  in  p-6)  and  (3-  #)  are  taken  over  *  semi-taEnlfe  frterval  bett-use  the 

99  m 

d^  te)  are  assumed  to  he  Awlytlc.  ’ 


Functions  -vith  tte  properties  (3-5)  aim  (3-6)  c; 
Taylor  coefficients  of  a  given  signal,  The  method  is  simply 
(3-1),  so  that 


n  he  M«d  ted 


hsHfefe  --5I  ir.  t  o 


F(fe) 


C  fa)  0  t&?)  /  d  (a 
■  •  n  *  0 


Both  sIms  of  (3-8)  can  he  nmBHolfed  far  C  to)  "  0  (-  d  #«> ,  • 

0  in 
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3.  l  — Continued, 


C  {mf 1  [F(cu)  0  (o>)  0*  (w)  ] 
o  m 


EC  (a)  f 

C(m)  n 
n=0 


C  0n  M  0^  M  3  • 

(3-9) 


Taking  the  Fourier  transform  of  both  sides  and  utilizing  property  (3-6), 


C(m) 

2a- 


F  (o>)  0q  (w)  (w)  e  du>  ss  f 


K  ^ 
m  0 


(t)  (3-10) 


m 


where 


j  0m<w>  1 2  eia,tdw 


(3-11) 


is  the  autocorrelation  function  of  the  time  signal  <b  (t)  whose  Fourier  transform  is 

mv 

0  (uO-  Since 
m 


max  R,  (t) 
t  m 


<*>  <w) 

m 


,  2 
I 


d  w 


§ 


1. 


we  have 


max 

4- 


F(o>) 


0  («)  0 
o  '  m 


(w) 


jwt  A 
e  qw 


J 


J 


fa 


m 


(3-12) 
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3. 1 _ — Continued. 

The  approximation  in  (3-12)  depends  upon  the  extent  to  which  the  inequality  (3-6)  holds 
true.  It  will  be  shown  in  Section  4,  however,  that  the  f  n  's  can  be  accurately  deter¬ 
mined  even  if  the  ’’much  less  than”  sign  in  (3-6)  is  replaced  by  a  simple  ’’less  than" 
sign. 


A  set  of  functions  that  satisfy  (3-5)  and  (3-6)  has  already  been  derived. 
The  signals  are  given  by 


*  (o>)  =  U(w/k")  /  E1/£  (3-13) 

n  n 

where  U  (co)  is  given  by  Equation  (2-8)  and  E  is  the  energy  of  the  function  U  (w/kn)  . 
Writing  the  left  hand  side  of  (3-9)  in  terms  of  U  (w). 


.  "1 


C  to)  F  (co)  U  (w)  U  (u  /  k  )  /(EE) 

o  m 


1/2 


_  C  (m)"1  F  (tv) 


<En  E  J 

o  m 


[u  <w>  c  (m)_1  w®  U  (co)  ; 


|C(m)|~  J™  [  i»m/'  2  U  <w)  J  [»nw'2U(»)]‘ 


(E  E  ) 
'  o  m 


1/2 


.-2 


LCJm)J__^_(co)  (  C(m/2)  {2  j  U(M/km/^  , 


(E  E  ) 
o  m 


(3-14) 


3-6 


N 


The  filter  transfer  functions  ]  U  (a?  /  km^  ")  { 2  ,  m  =  0,  1,  2,  , , . . , 
constitute  a  bank  of  constant  -Q  filters,  since  each  transfer  function  is  a  stretched 
version  of  the  preceding  transfer  function.  More  exactly,  the  ratio  of  center  frequency, 
a)  ,  to  centralized  mean  square  bandwidth,  D  -  o  ,  is  constant  The  constant-Q  prop- 

O  w  o 

erty  is  easily  demonstrated  by  defining 


=  / 

o  o 


O)  I  U  (w)  J 2  d  w  /  f  !  u  (w)  j z  d  u 

■Vi 


D  = 

U) 


»  « 
j  u2  J  U  (w)  j 2  d  u  /  f 


|  U(w)  |2dw]l/2,  (3-15) 


jj*  /  2 

3y  changing  variables  in  the  integrals,  it  iseasily  show's  that  when  U  {u>)— U  {&»/  k  *  ), 


a) 


knv'  2 

lv  W 


D 


km/2 


D 

0) 


and  Q  =  u)  /  (D  -  cs>  )  remains  unchanged. 
O  it)  o 
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3. 1  —Continued* 


The  above  mathematical  manipulations  can  be  summarized  as  follows: 

In  order  to  measure  the  spectral  coefficient  f  of  a  given  time  function 

f  (t),  process  f  {t>  with  a  filter  whose  transfer  function  is  J  U  (a?  /  km>/  )  |  ,  where 

U  (u>)  is  given  by  Equation  (2-8).  The  maximum  filter  response  is  proportional  to 

f  ,  provided  the  parameter  k  is  sufficiently  large. 
m 

To  measure  all  the  coefficients  f  ,  f  ,  ,  f  ,  use  a  bank  of  filters 

with  transfer  functions  j  U  /  k11^  )  j  ,  n=0,  1,  2,  ,  N.  This  bank  of 

filters  has  the  property  that  all  the  transfer  functions  have  the  same  Q,  i.  e. ,  the 
same  center  frequency  to  bandwidth  ratio. 

If  one  is  interested  in  the  spectral  Taylor  coefficients  of  a  filter  with 
transfer  function  F  (w)  *  one  can  pass  the  signal  u  (t),  with  Fourier  transform  U  (&>)  , 
through  the  filter.  The  resulting  signal  (with  spectrum  F  (w)  U  (a»))can  then  be  pro¬ 
cessed  with  a  bank  of  filters  whose  transfer  functions  are  proportional  to  U*  (cc  /  kn) , 
n  =  0,  1,  2,  ...  ,  N.  This  filter  bank  also  exhibits  the  constant  Q  property. 

If  one  were  to  find  a  constant-Q  filter  bank  in  the  laboratory  (or  In  nature), 
and  if  the  filter  transfer  functions  had  the  correct  shape,  one  could  conclude  that  the 
filter  bank  is  designed  for  Taylor  series  spectral  analysis.  The  filters  could  be  used 
for  characterizing  linear  systems  if  one  could  generate  the  proper  signal,  u  (t)  , 
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The  Ear  as  a  Spectrum  Analyzer, 


24  25  26 

As  a  con  sequence  of  masking  experiments  '  *  and  an  analysis  of 

27 

cochlear  dynamics  ,  the  human  ear  is  commonly  modeled  as  a  spectrum  analyzer 

composed  of  const ant-Q  filters.  The  exact  shape  of  these  filters  seems  to  be  dependent 

78 

upon  the  level  of  the  signals  that  are  used  in  the  masking  experiments.  For  reasonably 

low  levels,  however,  the  shapes  of  the  ear’s  filter  transfer  functions  are  very  similar 

2 

to  the  shape  of  j  U  (o>)  j  ,  where  U  (cj)  is  given  by  Equation  (2-8),  provided  f  G  (&?)  |  =1, 
That  is. 


U(w) 


jj  g-[(log  a.’)2  /  2  log  kj^  j  2  "  n  log  tii/  log  k 


(3-16) 


where  n  =  0,  ±1,  *2,  ...  ,  Consider,  for  example,  Zwicker’s  excitation  curve  fora 

1  kHz  signal  at  40  dB.  Zwicker’e  curve,  plotted  in  terms  of  voltage  level  (rather  than 

dB)  vs  frequency  (rather  than  critical  band  units)*?  is  shown  la  Figure  3-1,  Figure  3-2 

2 

is  a  scaled  version  of  the  function  exp  [  -32  (log  u )  J  .  The  two  curves  are  very  nearly 
identical. 


As  a  consequence  of  this  similarity,  it  can  be  postulated  that  the  ear  analy¬ 
ses  signals  In  terms  of  the  Taylor  coefficients  of  their  Fourier  transforms.  This  obser¬ 
vation  may  be  very  important  for  the  efficient  analysis ,  synthesis,  and  transmission  of 
epeech. 


There  is  one  disconcerting  aspect  concerning  the  similarity  between  the 


curves  in  Figures  3-1  aad  -3-2.  tt 


|  U  (a?)  | 2  =  exp  C  -32  (log  w)2 


(3-17) 


Figure  3-1,  Zwtoker'a  Excitation  Curve  for  a  Signal  Frequency  of  1  KHz  at  an  Excitation  Level  of 

40  da  The  dfmpe  of  this  Curve  Is  Similar  to  Masking  Data  Measured  by  other  Research- 
era.  24.  2,>  Note  Linear  Level  Scale  (rather  than  dB)  and  KHz  frequency  scale  (rather  than 
critical,  tod  units). 


Figure  3-2.  A  Plot  of  j  U  <cu)  p  •  oxp  C  -32  {log  w)M 
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3,2  —Continued. 


then  log  k  =  1/32  .  or 


k  =  1,0318  .  (3-18) 

When  $  (&?)  is  defined  as  in  (3-13),  the  left  hand  side  of  (3-6)  becomes 

m  _  tt  g  £ 

max  1  x,JU  (t,  k  )  |  ,  where  {  XIJU  (T,  s)  |  is  the  wideband  ambiguity  function 

of  the  signal  with  Fourier  transform  (2-8),  It  can  be  shown  (Appendix  and  Reference  3) 


that 


m3X  ,  *uu  (t*  ^m~n)  I2  =  *W/2  . 

t  UJ 


(3-19) 


For(m  -  n)=  1,  Equation  (3-6)  can  only  be  written  with  a  simple  "less  than"  sign  (rather 
than  << )  for  the  value  of  k  given  by  (3-18).  This  situation  is  not  disastrous,  however, 
because  a  linear  transformation  of  the  filter  outputs  can  still  result  In  determination  of 
the  individual  coefficients  f  m.  Such  a  transformation  arises  naturally  when  a  likelihood 
ratio  test  is  applied  to  the  filter  outputs,  A  receiver  that  implements  a  likelihood  ratio 
test  is  discussed  in  Section  4, 


It  should  also  be  mentioned  that  the  bandwidths  of  the  ear’s  filter  bank  may 

29 

be  increased,  allowing  for  a  larger  value  of  k.  *  D.  M,  Green  has  found  that  human 
listeners  increase  their  effective  bandwidth  to  detect  a  wideband  noise- like  signal  masked 
by  additive  rsoiae.  Green’s  data  suggests  that  "critical  bands  are  not  fixed  in  width  but 
are  adjustable  so  as  to  match  the  particular  detection  situation. " 


*  It  is  shown  in  the  Appendix  that  the  centralised  bandwidth  Is  proportional  to  the  square  root  of 
k3/  2  (k  1/  2  -  1),  A  larger  value  of  k  thus  tellies  a  larger  centralized  bandwidth,  and 
conversely. 
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3, 3  Animal  Echolocation, 

The  foregoing  results  provide  an  extremely  simple  and  elegant  approach 
to  our  previous  observations  about  animal  echoloeation  signals  (Section  2),  Suppose  that 
a  given  receiver  analyzes  auditory  signals  in  terms  of  the  Taylor  series  coefficients  of 
their  Fourier  transforms.  If  such  a  receiver  incorporates  a  constant  Q  (or  critical 
bandwidth)  filter  bank,  it  effectively  performs  the  operation  shown  in  Equation  3-14. 
According  to  the  left  hand  side  of  3-14,  incoming  data  is  multiplied  by  U  {&?),  where 
U  (cu)  Is  given  by  (2-8).  The  resulting  product  is  then  multiplied  by  U*  (a?  /k*  )  , 
m  =  0,  1,  2,  ...  . 


Consider  the  utilization  of  such  a  receiver  in  a  sonar  system.  As  before, 
it  is  assumed  that  targets  can  be  characterized  In  terms  of  their  transfer  functions. 
Letting  F  (w)  represent  the  target  transfer  function,  a  target  can  be  characterized  by 
transmitting  a  signal  U  (&>),  as  in  Equation  (2-8).  The  echo  is  then  the  product  U  (as)  X 
F  <co).  This  product  should  be  multiplied  by  U  *  (u/k*”)  ,  m  =  0,  1,  %  . . . ,  as  in 
Figure  1-1,  in  order  to  determine  the  power  series  coefficients  o  F  fo). 


Given  a  set  of  constant-Q  (critical  bandwidth)  filte.s.  a  natural  method  of 
signal  analysis  is  to  determine  the  Taylor  coefficients  of  received  signal  spectra.  Sonar 
target  description  can  then  be  accomplished  by  determining  the  Taylor  coefficients  of 
target  transfer  functions.  The  waveforms  that  should  be  used  for  such  a  characteri¬ 
zation  are  the  set  of  signals  that  have  already  been  matched  to  those  employed  by  the 
Atlantic  Bottlenose  Dolphin  l  Filtered  versions  of  the  same  signals  are  used  by  the 
bat  Myotis  lucifagus.  TMs  observation  is  further  <Hscussed  in  Section  6. 
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3. 4  Discussion, 


Our  ntodei  of  ear  as  a  particular  kind  of  spectrum  analyzer  justifies  tte 
use  of  a  specific  class  of  functions  for  animal  echolocation.  The  fact  that  these  func¬ 
tions  are  actually  utilized  bv  oats  and  dolphins  aids  further  validity  to  our  model  of  the 
auditory  processor. 


The  above  analysis  strongly  suggests  that  the  class  of  functions  (3-16)  is 
suitable  not  only  for  utilization  by  animals,  bat  by  humans  as  well.  Indeed,  when  a 
dolphin- like  pulse  is  reflected  from  different  targets  and  suitably  scaled,  people  can 
tear  the  difference  between  the  echoes.  This  experiment,  together  with  the  fore¬ 
going  analysis,  suggests  that  the  signals  should  be  incorporated  into  ultrasonic  guidance 
devices  for  blind  people"0  and  into  small  &.  s  for  scuba  divers. 31 


Finally,  the  above  results  suggest  that  a  study  of  animal  echolocation  can 
provide  important  clues  for  analysis  of  speech,  both  of  cetaceans  (if  their  communication 

sounds  can  indeed  be  classified  as  speech)  and  unmans.  These  ideas  will  be  reconsidered 
In  Section  1 1. 


ENT  A 


;i.IHOOD  RATIO  TESTING  WITH  FILTERS  THAT  1MPLEM 
TAYLOR  SERIES  SPECTRUM  ANALYSIS. 


4. 1 _  Processing  the  Filter  Responses. 

The  receiver  in  Figure  1-1  does  not  felly  specify  the  processing  that  should 
be  applied  to  the  outputs  of  the  constant  Q  filters.  A  more  complete  specification  can  be 
obtained  by  subjecting  the  filter  cuQmts  to  a  likelihood  ratio  test. 

The  filter  outputs  are  to  be  processed  in  order  (i)  to  detect  a  signal  with 
known  Fourier  transform  F  (tu)  or  (111  to  estimate  the  Taylor  coefficients  of  F  (tc)  when 
they  are  not  known  a  prion,  i.  e. ,  to  perform  Taylor  series  spectrum  analysis. 


When  the  input  signal  has  Fourier  transform  F  (w)  U  (&?),  it  will  be  shown 
that  the  same  initial  operation  is  performed  on  the  receiver  outputs  ft)  to  decorrelate 
the  noise  responses  for  maximum  likelihood  detection  of  a  known  spectrum  F  fee),  and 
(i  Ij  to  determine  the  Taylor  coefficients  of  an  unknown  spectrum  F  (art.  It  will  be  shown 
that  the  ideal  detector  implements  a  simple  correlation  process  on  the  filler  outputs, 
while  the  maximum  likelihood  estimator  measures  the  mean  of  a  transformed  version 
of  the  filter  responses. 

In  the  case  of  sonar,  U  {&?)  is  the  transmitted  signal  spectrum,  F  (cc)  Is 
the  target  transfer  function,  and  U  (w)  F  (w|  is  the  echo.  When  F  {&>)  itself  Is  the  re¬ 
ceived  signal,  one  can  determine  the  Taylor  coefficients  of  the  spectrum  H  |w),  where 
F  C&i)  =  H  («)  U  (<■>).  Since  the  function  U  (w)  has  known  coefficients,  the  coefficients 
of  F  (bJ)  can  be  determined  from  the  measured  H  (w)  coefficients.  R  mill  be  shown  that 
the  F  (u>)  coefficients  a*e  obtained  from  those  of  H  (to)  by  means  of  a  convolution  operatic 


s  with 


Assuming  that  the  received  signs!  has  a  spectrum  F  {fa?)  U  (fa?)  or  F  (w)  s 
H  {»)  U  (&?),  we  want  to  determine  the  Taylor  coefficients  of  F  (fa?)  or  H  {fa?}.  According 
to  {3-14),  the  signal  can  be  analyzed  with  a  bank  of  filters  with  transfer  functions 


U*  {«  /  k‘“l  /  C  tot  (EE  i*7  nt=0,  1,  2,  ...  .  To  determine  the  weighted  coef- 

O  13 

fieients  c  (m)  f  ,  the  fifier  bank  consists  of  unit  energy  transfer  functions 
m 


Z  (fa?)  =  U*  {fa?  /  km)  /  E  2/  2 
m  m 


/4-1) 


where  E  si. 
o 


N  <t> 
in 


noise  at  output  of  filter  number  m 
impulse  response  of  filter  number  m 


inverse  Fourier  transform  of  Z  (fa?) 

13 


N  (t> 

<l> 


noise  at  the  input 


autocorrelation  function  of  N  (t). 


The  likelihood  ratio  test  involves  a  computation  of  two  coi-arinr.ee  matrices 


covariance  matrix  at  the  filter  outputs,  one  most  determine  E  j  N  *)  N  *  (t)  l 

*  n  tt=  * 
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4. 3 _ —Continued. 

Letting 


*  I/O 

f  =  f  C  (n)  E  '  , 

n  n  n 


and  using  the  relation  {see  Appendix) 


(4-8) 


*uu  <°’  k 


n-ni)  _  -(n-  m)  /4 


(4-9) 


Equation  (4-7)  gives 


g. 


g* 


-1/4 


-1/4 


k-NZ/4  k-(N-l)2/4 


*  «  •  • 


,-N2/4 


-{N-l)2/4 


k' 


J  • 
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4. 3 _ — Continued. 


Letting  be  the  covariance  matrix  of  the  noise  outputs,  Equations  (4-5),  (4-6),  and 
(4-9)  imply  that 


K 

=o 


E  (N 


N 


m 


N 


-(n- 


m)2  / 4 


(4-11) 


or 


g 


(4-12) 


Another  covariance  matrix  must  be  computed  in  order  to  apply  a  likelihood  ratio  test. 

This  matrix  is  the  covariance  of  signal  plus  nGise  at  the  filter  outputs.  For  noise  with 
zero  mean,  the  elements  of  the  matrix  are 

E  {  K+n„<'>]  [■£♦*;«]}  -  E{*„}  E{g-} 

*  E  {  Nn  (t)  N4  '-U> 

where  it  has  been  assumed  that  the  components  g  correspond  to  a  deterministic  signal, 
sothat  E  {gngm*}  =  E  (g„)  E  jgm*|  . 

Letting  5  be  the  covariance  matrix  of  'liter  outputs  when  a  nonrandom 
signal  and  additive  noise  are  present  at  the  input  of  the  receiver.  Equations  (4-13)  and 
(4-11)  give 

*In  Equation  (4-11),  the  symbol  C  f  (a,  m)  3  indicates  a  matrix  such  that  the  element  in  the 
ntb  row  and  mtk  column  is  given  by  f  (n,  m). 
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4.3  — Continued, 


k-(n-m)“/4 


(4-14) 


Assuming  that  tbe  input  noise  N  <t>  is  Gaussian,  the  likelihood  ratio  test 
compares  the  quantity 

l  (B)  =  (BT-  M  T)  Q  (R-  M  )  -  (RT  -  M  T)Q  (R-  M  )  (4-15) 

0=0  O  1-1  1 

with  a  threshold.  Si  (4-15),  R  is  the  data  vector,  i.  e. ,  the  voltages  that  appear  at  the 
filter  outputs  at  any  given  time.  The  vector  M-  is  the  mean  of  the  data  vector  when  the 
signal  is  absent,  and  M.  is  the  mean  of  the  data  when  signal  is  present.  For  the  case 

*  -X  _1 

under  consideration,  M  =0  and  M,  =  e .  Finally,  Q  =  K  and  Q,  =  K„  , 

o  1  — o  =o  —1  =1 

For  the  case  of  a  nonrandom  signal  in  the  presence  of  white,  Gaussian 

noise,  we  have  determined  that  K  =  K  .  in  accordance  with  Equation  (4- 14).  Li  such 

32  u 

a  case,  (4-15)  simplifies  to  the  expression 

I  (B)  =  B  T  Q  <M,  -  M  )  =  RT  Q  g.  .  (4-16) 

— O  a  o  — o 


The  above  expression  is  very  meaningful  in  terms  of  Equations  (4-10)  and 

(4-11),  since 


f 


N 

o 


Q  g 

= o  — 


(4-17) 
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4.3  — Continued. 


Substituting  (4-17)  into  (4-16)  and  neglecting  the  constant  multiple  N^/  2  (since  a  constant 
multiple  merely  changes  the  value  of  the  threshold), 

i(B)  *  5  i  =  y  r  f  .  (4-18) 

jL  °  n 
n=0 

Equation  (4-18)  describes  a  correlation  process.  The  filter  outputs  r 

A  tl 

are  correlated  with  the  signal  coefficients  f  .  The  process  is  illustrated  in  Figure  4-1. 

n 

In  (4-16),  the  receiver  operates  on  the  data  vector  R  In  such  a  way  as  to 

decorrelate  the  noise  readings.  The  transformed  data  should  then  be  correlated  with 

the  expected  values  of  the  filter  outputs  when  the  signal  is  present.  These  values,  tow- 

1/  2 

ever,  are  related  to  the  weighted  coefficients  C  (n)  f^  by  a  transformation  which 

is  the  inverse  of  the  first  (decorrelation)  transformation.  The  result  is  a  direct  cor¬ 
relation  of  filter  outputs  and  known  coefficients,  with  no  matrix  transformations  at  alL 


4. 4 _ The  Effect  of  Linear  Filtering  on  the  Taylor  Coefficients  of  a  Signal's 

Fourier  Transform. 

It  has  been  assumed  that  (i)  the  input  noise  is  white  (Equation  4-4)  and 
(I  I)  that  the  desired  coefficients  are  those  of  F  (w),  whew  the  input  signal  has  spectrum 
F  (sj)  U  (u>).  With  these  assumptions,  a  simple  processing  implementation  has  been 
obtained.  In  the  interest  of  simplicity,  it  is  therefore  reasonable  to  assume  (i)  that  a 


4-8 


4-1,  Detection  of  a  Signal  with  Spectrum  F  (w)  U  (to)  When  the  Taylor  Series  Coefficients 

of  F  (w)  are  known. 
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4.4  — Continued. 


whitening  filter  precedes  the  receiver  when  colored  noise  is  present  and  (i  i)  that  all 
input  spectra  can  be  written  as  the  product  of  U  {<*?)  with  a  finite  order  polynomial 


n=0 


f  « 
n 


To  discover  the  effect  of  these  assumptions  on  the  measured  coefficients, 
it  is  necessary  to  investigate  the  change  in  a  spectrum's  Taylor  series  when  the  corre¬ 
sponding  time  signal  is  passed  through  a  linear  filter. 

Consider  the  product  A  (a?)  B  (u),  where  A  (o>)  can  be  visualized  as  the 
Fourier  transform  of  a  signal  which  is  processed  by  a  linear  filter  with  transfer  func¬ 
tion  B  (a?).  Let 


A  (a?) 


2 

a  +  a,  fa?  +  te  + 
o  1  2 


a..  fa> 
?» 


(4-19) 


and 


2 

B  (fa?)  =  b  +b  u  +  b  w  4 . 

O  1  2 


<4-20) 


*  Notice  that  y  f  fa?”  is  neither  tend  limited  nor  energy  limited  (square  integrable) 
n 

n=0 

for  N  bounded  and  0  s  a?  s  «  ,  For  U  (a?)  given  by  (2-8),  however,  tte  function 

n  -  2 

T  j  f  a?  3  U  (fa?)  is  both  tend  limited  (in  tte  mean  square  sense)  and  in  L  (0,  * ). 

'  n-0  n 


I 

I 


0 

8 

0 
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4.4  — Continued. 


Now  consider  the  convolution  of  two  signals  a  (s)  and  b  (x)  ; 


c  (r ) 


a  {r  -  x)  b  (x)  dx  s  a  (x)  *  b  (x)  . 


(4-23) 


Tha  convolution  of  the  vectors  (a  ,  a,,  . .  *  a  i  mid  (b  ,  b,  ....  bj  is  then  obtained  by 

o  1  IT  o  1  N 

reversing  the  a- vector  and  translating  it  with  respect  to  the  b-  vector.  For  each  trans¬ 
lation,  the  convolution  is  the  sum  of  the  products  of  the  corresponding  elements.  In 
matrix  notation. 
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spectra. 


Since  the  coefficients  of  the  product  of  two  spectra  can  be  obtained  by  a 
convolution  operation,  they  can  also  be  obtained  by  a  multiplication  of  Fourier  trans- 
forma.  One  can  construct  pulse  trains  of  the  form 


a  (x)  =  \  a  reel  (x-  2n) 
P  Zf  n 

n=G 


b  (x)  =  V  b  reel  (x-  2n) 

P  L  n 


{4-25) 


where 


rect  (x)  = 


0  £  x  -r.  1 


0  ,  x  <  0,  x  >  1 


%  (x>  =  %  (*>  *  M*) 


c„  tri  (x  -  n) 
u 


y;i 


4, 4  —  Continued, 


where 


iri  <x) 


x  ,  Osssi 

2  -  x,  lsx^2 

0  ,  x  <  0,  x  >  2  . 


{.1-28) 


The  coefficients  e  in  (4-27)  are  given  by  (4-24),  To  obtain  these  coef- 
n 

flcients  more  easily,  one  can  compute  Ap  {&?)  and  Bp  (tu) ,  the  Fourier  transforms  of 
ap  (x)  and  bp  (x) .  The  function  Cp  (x)  is  the  inverse  Fourier  transform  of  Ap  (a?)  x 
Bp  (a*). 


As  in  the  more  conventional  analysis  of  linear  systems  utilizing  sine  waves, 
Fourier  transformation  Is  seen  to  be  a  valuable  toot  Usually,  Fourier  transforms  are 
multiplied  in  order  to  compute  convolutions  that  describe  time  functions.  In  the  above 
analysis,  Fourier  transforms  are  again  utilized  to  compute  convolutions.  For  Taylor 
series  spectral  analysis,  however,  the  convolutions  describe  functions  in  the  frequency 
domain.  Convolution  of  two  time  functions  is  equivalent  to  convolution  of  their  Taylor 
series  spectral  coefficients. 


4. 5  Detection  of  a  Signal  F  (u>)  with  Known  Coefficients. 


In  order  to  detect  a  known  signal  F  (w)  (rather  than  F  (w)  U  (£c»wtth  spectra 


coefficients  £,  Equation  (4-18)  implies  that  the-  filter  responses  x.  should 
with  a.  i«o.  » 


be  correlated 


n^O 


11  n 


where,  from  (4-  8), 


b 

n 


C  (n) 


1/2 


(4-30) 


so  that 


MS  =  £  rn  CHE/’  hn  . 
n-"0 

The  vector  h  corresponds  to  the  coefficients  of  H  {«),  where 


(4-31) 


F  {ii)  s  H  («)  U  (w)  , 


If  U  (0)  has  coefficients  ti .  then  according  to  (4-2S). 


inimMi 


I 

I 
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4, 5  — Continued, 
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(4-33) 


uhere  the  approximation  Is  obtained  by  representing  U  (ot)  in  terms  of  a  polynomial  with 
N  coefficients  on  a  bounce d  frequency  Interval. 


Notice  that  the  form  {4-33}  would  net  exist  If  U  (w)  were  defined  In  terms 
of  sn  infinite  number  of  coefficients.  When  U  (ca)  is  approximated  by  a  finite  number  N 
of  coefficients,  the  frequency  interval  over  which  the  approximation  Is  accurate  must 
always  be  bounded.  ♦  3y  the  same  argument,  if  the  input  function  F  (ce)  Is  square  fate- 
grable  and  is  deocrioed  In  terms  cf  N  coeffitteitts  «W  bounded),  the  analysis  must  apply 
te  a  ecu  need  frequency  Interval  ratter  than  to  the  semi-infinite  interval. 


W 
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4, 5 _ — Continued, 

Equation  (4-33)  can  be  written 

l  ^  a  *  ii  =  'a,  *  (4-34) 

where  T  is  the  square  (SxN)  matrix  in  (4-33).  The  filter  outptfs  in  Equation  (4-31) 

1/5 

should  therefore  be  correlated  with  C  (n  l  E  "  h  ,  where 

n  n 

Jl  &  T  _1  1  (4-35) 

=u 

and  f  corresponds  to  the  known  coefficients  of  the  input  signal  spectrum  (or  approxima¬ 
tions  to  these  coefficients  that  are  valid  over  a  bounded  frequency  interval). 

The  process  Is  shown  to  Figure  4-2.  In  terms  m  Figure  4-1,  F  (ai)  is  writ¬ 
ten  as  H  {«)  U  (tit).  The  receiver  then  works  to  detect  H  (sy)  as  in  Figure  4-1,  where 
H  (u)  =  f  u  (a)  T1  F  (w)  or  h  =  T  _1  X- 


4.6  _  Estimation  of  Unknown  Spectral  Coefficients. 

If  the  signal  coefficients  are  unknown,  they  can  be  estimated  by  determla- 


!i<jiir.i||i|l||.pii  iiiriiM.ii|iiiiiii|i||ii||i|.iiiiiii|i| 
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4.6  —Continued. 


describing  the  filter  outputs  in  the  presence  of  a  signal,  the  coefficients  g  (Equation 

n 

4-7)  can  be  estimated  by  setting 

loSCPfSI  H.)  J  -  0  .  (4-36) 

8gn  1 

33 

Equation  (4-36)  leads  to  a  maximum  likelihood  (ML)  estimate  of  the  g  coefficients. 

Once  again,  the  fact  that  f. is  determined  from  £by  computing  Qq  ^(Equa¬ 
tions  4-17  and  4-8)  results  in  considerable  simplification.  Consider  an  ML  estimate 
performed  on  Qq  r,  rather  than  on  r  itself.  The  transformation  Qq  causes  the  noise 

readings  to  be  statistically  independent.  In  the  presence  of  signal,  Qq  B.is  a  vector  of 

** 

independent  Gaussian  random  variables  with  means  given  bv  Q  g  =  (2/  N  )  f  ,  from 

=  o  —  o  — 

Equations  (4-7)  and  (4-17).  The  distribution  p  (Qq  B.  |  H^)  is  therefore  maximized  by 
determining  the  mean  value  of  each  element  of  the  output  vector  r  ; 

E  {90  51  Hjl  =  (2/No)  £ 

- ►  M  L  estimate  of  f  ,  (4-37) 


32 

*  The  notation  used  in  this  section  is  from  Van  Trees'  book  .  By  using  capital  R  and 
lower  case  r.,  Van  Trees  dlstinquishes  between  a  random  variable  (B)  and  sample 
values  of  a  random  process  (i). 
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4,0  — Continued. 


where,  from  (4*  fy , 


f  =  [  C(n)  E  1/2  f1  f  .  (4-38) 

n  n  n 

The  above  process  is  illustrated  in  Figure  4-3. 

The  processor  in  Figure  4-3  estimates  the  coefficients  of  F  (w),  where 
the  received  signal  is  F  (w)  U  (w).  If  the  received  signrt  is  F  (a?)  alone,  one  can  use 
the  same  procedure  to  estimate  the  coefficients  of  H  (u>),  where  F  (e>)  -  H  (w)  U  \o})  . 
To  obtain  the  coefficients  of  F  (cj)  from  those  of  H  (a>),  Equation  (4-34)  can  be  invoked. 
The  resulting  process  is  shown  in  Figure  4-4  . 


4, 7  Relaxing  the  Decorrelation  Requirement. 


Equation  (4-17)  implies  that  the  Taylor  coefficients  of  an  incoming  spectrum 

can  be  estimated  if  the  K  matrix  (the  covariance  matrix  of  the  noise  responses)  can  be 

-o 

inverted.  Inversion  of  K  is  possible  if  its  determinant  is  nonzero,  i.  e, ,  if  the  para¬ 
mo 

meter  k  in  (2-3)  and  (4-10)  is  greater  than  one.  In  other  words, 

k<n-m)2/4  .  Xuu  (  0»  kn  ”  m  )  <  1  .  (4-39) 

It  has  already  been  remarked  that  complete  decorrelation  of  filter  transfer 
functions  is  unnecessary.  It  is  not  necessarily  required  that  [  xm  (0  ,  kn  m  )  |  "  be 
much  less  than  one  ;  a  simple  "less  than1'  sign  is  sufficient. 


Figure  4-4.  Maximum  Likelihood  Estimation  of  the  Taylor  Series  Coefficients  of  F  (w),  when 

the  Input  Signal  is  F  (to). 


ESL-PR115 


4.8 


Detection  of  a  Signal  with  Spectrum  F  (a?)  U  (up  when  the  Power  Series 
Coefficients  of  F  (a?)  are  Random  Variables.  ~ 


The  receiver  in  Figure  4-1  is  easily  generalized  to  the  case  where  each 

coefficient  f  is  a  random  variable  with  its  own  probability  density  function  (pdf).  Target 
n 

coefficients  that  are  random  variables  can  result  from  the  aspect  dependence  of  a  given 
reflector  and/  or  the  differences  between  in  lividiial  members  of  a  single  target  set,  e.  g. , 
the  set  of  dangerous  sharks. 


Unless  a  reflector  is  extreme  y  symmetrical,  its  transfer  function  F  (a?) 

will  be  dependent  upon  the  position  from  which  the  reflector  is  observed.  The  Taylor 

series  coefficients  of  F  (cj)  are  therefore  dependent  upon  aspect  angle  ,  6.  Let  F{fJ0) 

be  the  pdf  describing  the  aspect  dependence  of  each  coefficient  For  a  particular 

target,  P  (f  i  9)  equals  S  f  f  -  f  (8)  ],  where  f  (0)  is  a  deterministic  function  deserib- 
n  n  n  n 

ing  the  aspect  dependence  of  f  .  For  a  set  of  similar  but  not  identical  targets,  P  (f  10) 

n  n’ 

has  a  nonzero  variance.  Let  P  (0)  be  the  probability  of  encountering  the  target  from  a 

particular  aspect  angle,  0 .  Then  the  pdf  associated  with  the  random  coefficient  f  is 

n 


P,y  = 


£i 

j  P(fn  !  8)  P  (6)  d  9 


(4-40) 


Figure  4-1  summarizes  the  analysis  In  Section  *1-3,  where  a  likelihood 

ratio  test  is  applied  to  a  receiver  with  filters  Z  (w)  =  U  (w  /  k  )  /  E  1  "  .  The 

m  1/  ^ 

test  consists  of  correlating  the  filter  outputs  with  the  numbers  |c(n)  E^  "f  |  ,  where 
If  |  is  a  set  of  H  power  series  coefficients  describing  F 


♦The  aspect  angle  8  is  a  scalar  if  the  problem  is  considered  in  two  dimensions;  8  is  a  two 
term  vector  for  three  dimensional  problems. 
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4. 8  — Continued, 


When  the  coefficients  J  f  j  are  random  variables, 


/ 


Hm  =  I  £{R  i  D  P<0  df 


where,  from  (4-18). 


Therefore, 


£  (B  i  D  = 


N 


n=0 


1/2 

C  (n)  E  J/  r  f 
n  n  n 


where 


N 


£(B)  « 


■  2 
n®0 


C(n)  E  1/2  r  f 
n  n  n 


f  = 
□ 


a 

/ 


f  P(f)  df 
n  n  n 


Substituting  (4-40)  into  (4-44)  , 


7° 

it  can  be  shown  “  that 

(4-41) 

(4-42) 

(4-43) 

(4-44) 
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--Continued. 


®  2  ~ 

1 1 


f  P  (f  !  6)  P  (6)  d  &  df 
n  n  ‘  n 


-O  O 


2  » 

j  s  jfn  {  ej  P(8)dP 


<4-45) 


where  the  conditional  expected  value  is 


E  !fn  1  e\ 


/  fnp<f„i(’>  d,„  • 


(4-46) 


If  the  problem  is  to  detect  a  single,  well- specified  target,  E  jf  |  0 }  is 
just  the  function  (6)  describing  the  dependence  of  each  coefficient  upon  aspect  angle. 
In  this  case, 


f  (6)  P  (0)  de 

3 


1  C»  EnI/2  /'  fn 


(P)  P  {©)  d  P 


WM . .  am 
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5. _ SIGNAL  TO  INTERFERENCE  RATIO  (SIR)  MAXIMIZATION  USING 

LINEAR  DETECTION, 


5, 1 _  Linear  Defection  of  Received  Signals. 


Is  radar  systems,  decisions  are  usually  based  upon  an  envelope  detected 

or  squared  v:  rsion  of  correlator  response.  The  reason  for  such  nonlinear  processing 

?6 

is  the  elimination  of  "fine  structure"  within  the  autocorrelation  function.  This  fine 

structure  results  from  modulation  of  the  baseband  signal.  If  the  signal  u  {t}  modulates 

a  carrier  with  frequency  u  ,  the  autocorrelation  function  of  the  modulated  signal  has 

o 

the  same  envelope  as  the  autocorrelation  function  of  u  ft).  The  phase  of  the  autocorrela¬ 
tion  function  varies  as  w  t .  This  phase  variation  is  the  so-called  "fine  structure"  which 

o 

is  eliminated  In  order  to  feed  a  "cleaner"  signal  into  the  threshold  detector. 


An  important  difference  between  wideband  sonar  and  radar  (or  narrowband 
sonar)  is  that  the  baseband  signal  Is  actually  transmitted  without  modulating  a  carrier. 

All  operations  can  therefore  be  carried  cart  at  baseband  frequencies,  and  the  autocorrela¬ 
tion  function  (for  a  properly  designed  signal)  is  no  longer  corrupted  by  the  "fine  structure 
of  radar  theory. 


Decisions  can  therefore  be  based  upon  the  response  of  a  completely  linear 
system.  Complete  linearity  is  advantageous  because  phase  differences  can  be  exploited 
for  both  clutter  suppression  and  enhancement  of  signal  to  noise  rs'io. 

For  clutter  suppression,  one  s.-in  design  a  filter  so  tost  tie  maximum 
response  to  an  unwanted  reflector  echo  is  negative,  whereas  the  maximum  response 
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5. 1  — Continued. 


to  the  desired  target  is  positive.  Square-law  detectors  would  require  the  unwanted  re¬ 
flector  to  invoke  a  filter  response  whose  magnitude  was  small. 

For  noise  suppression,  one  can  ignore  large-negative  Alter  outputs  if  the 
response  to  the  desired  signal  is  positive.  Half  of  the  noise-induced  peaks  of  the  squared 
filter  response  are  thus  ignored.  As  one  would  expect,  the  resulting  signal  to  noise 
ratio  is  doubled,  just  as  in  comparisons  of  coherent  vs,  noncoherent  detection. 

S.2 _ Signal  Design  for  Linear  Detection. 

If  linear  detection  is  to  be  used,  Ate  autocorrelation  fundi.—  should  ideally 
have  only  one  positive  peak.  It  would  seem,  however,  that  large  negative  sMelobes  are 
permissible.  The  large  negative  sidelcbes  of  the  autocorrelation  functions  *  in  Figures 
5-1  and  5-2  will  not  not  induce  false  triggering  or  range  aiebiquities  if  linear  detection 
is  used.  These  autocorrelation  functions  correspond  to  k  =  1.5  and  k  =  2  in  Equation 
(2-S)  ,  where  J  G  &)  J  =  1.  These  k-values  bound  most  of  the  animal  signals  that 
have  teen  studied  thus  far.  In  Tables  2-1  and  2-2  ,  for  example  ,  k  **  1.  G5. 

Although  large  negative  sidelohes  are  permissible  from  toe  viewpoint  of 
range  ambiguity  peaks,  such  sMelobes  are  apparently  not  desirable  for  the  resolution  of 
two  closely  spaced  targets.  A  small  positive  peak  can  be  masked  by  a  larger  negative 
peak,  corresponding  to  a  large  nearby  reflector.  The  peak  to  sidelcbe  ratio  improves 


O  (M) 


filtering  with  scaled  versions  of  the 


soeetmm,  as  mseasseti  ia  3 


When  linear  detection  is  used  In  contoarfiim  with  tise  signals  p-»),  ite 


one  complex  target.  The  existence  of  two  targets  is  then  determined  by  Taylor  serf 
spectral  analysis.  From  this  viewpoint,  signals  with  large  negative  2a;ac-.-rrel3.tica 
sidelobes  are  a  sain  remissible. 
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— Continued, 


and  where  the  parameter  k  is  sufficiently  large  30  that 


,  .  n  -  m  .2 

max  v  (T>  <  ) 

*uu 


nm 


(5-2) 


The  filter  is  modelled  as  a  bank  of  matched  filters  with  transfer  functions 


Z  (w) 
n 


w11  U  (co)  /  E 


(5-3) 


A  weighted  sum  of  the  filter  outputs  is  used  to  determine  the  presence  of  a  target, 
shown  in  Figure  5-4.  The  weights  v_  a  re  constrained  so  that 


N 


n=0 


i  v  | 
n 


2 


The  target  transfer  function  is  defined  as 


(5-4) 


F(w)  - 


N 

2  ^ 


11=0 


(5-5) 


and  the  average  clutter  transfer  ft  met  ion  is 


C  (a,*) 


The  target  echo  is 


c  u 
n 


n=0 


(5-6) 


D—  / 
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—  Continued. 


N  N 

F  (a*  U„.  <w)  =  V  Y  f  u  (J1  +  m  U  (a*) 

Trans  L*  Zi  R  m 

n-0  m=0 


f  u  U  lea)  *■  ({  u,  +  fu  )  u?  U  (a>) 

A  A  *  ft  J  7  n 


+  (fo  *2  +  fl  “l  +  f2  V  “  U 


The  clutter  echo  is 


N  N 


C  U rnu..'  <“’>  ■  2  2  '»  ■m"”"  U  <“> 

n=0  m=0 


=  c  u  U  (£J)  4  (c  U  fC«)  u->  U  (a?) 
oo  ox  lo 


+  (C  Uft  +  C  U  +  C  U  )  CU  TJ  {&))  -f  . 
O  a  11  4  0 


The  weighted  sum  of  the  filter  outputs  is 


(f  u  )  v  4  if  u,  *  f  u  )  v-  +  (f  +  f,  u.  +  f  u  }  v  +  . 
o  o'  o  %  o  1  1  o'  1  o  2  li  2  o  2 
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— Continued. 


for  the  target  and 


(c  u  )  v 
"  o  o  c 


<Vl 


j  u  )  V  i  <e  a  -r  c  a  f  c  u  t  v 
1  o  1  o  2  11  2  o '  2 


for  the  clutter. 


To  maximize  the  signal  to  clutter  ratio,  one  can  maximize  the  difference 
N  ,  2  2 

h  comraxnts  on  j  j  and  ^  1  v_  I  .  In  otter  words,  one 


-  Z  with  comraints  on 
1  c 


can  maximize  the  quantity 


(«»  2)  =  ZT-IC  -  M  j  I  un  12+  £  I  vn  |2)  .  (5-lx) 


Differentiating  J  (u  ,  y)  with  resect  to  v  and  setting  tte  result  equal  to  zero  gives 

n 


y  (f  .  -  c  ,)  u  -  2  A  V  *  =  0 
**  i  ii  n 


i,  j 

i  +  j  =  n 


2  X  Z 


■j  =  n 
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5. 3  —  Continued, 


In  other  words, 


u  -  A  v  *  and  v  =  A  u* 


where  A  is  the  NxN  matrix  in  (5-13).  It  follows  from  (5-15)  that 


T  * 

u  =  A  A  u 


and 

v  =  A  A  Tv  (5-17) 


f  »  *'T" 

so  that  u  and  v  are  eigenvectors  of  the  transformations  A*  A  and  A  A  respectively. 

For  example,  suppose  that  we  wish  to  detect  a  planar  reflector  with  F(«)  =  l 

g 

amid  clutter  with  C  (t*j)  -  co  .  Then,  for  N  -  5, 

10  0  0 

0  10  0 

,  0  0  10 
-10  0  1 


3 


3 


T  * 
A  A 


An  eigenvector  of  both  (5-19)  and  (5-20)  Is 


u  = 


Notice,  however,  that  the  s mm  matrices  A“  A  sr.d  A  A  would  result  If  F  (es)  -  cs™ 
a»w  C  -  1.  One  must  therefore  check  to  see  whether  the  SIR  has  been  maximized  oi 
minimized.  A  simple  check  is  to  eoraputeT^and 


I 


L-  P  Til  1 5 


of  a  Basis  Set  of  Wideband  Functions  over 


assume  that  toe  ranges  of  target  ana  clatter  are  knosr.  a  priori.  Alternatively,  one 
can  follow  the  usual  practice  of  envelope  defecting  the  filter  outputs. 


A  second  discouraging  feature  of  sinusoids  is  '.hat  they  are  sensitive  to 
iargr:  mottos.  Any  attempt  at  a  fine- grata  analysis  ©f  echo  spectr*  mast  to  long 
carattaa  siMBotds  which  are  serious;  .'  affected  fey  unforeseen  motion  of  target  or  chi- 
ter.  Hsfc  wideband  signals  |2-8i  are  generally  insensitive  to  small  stretch  factors  In¬ 
duced  by  reflector  motion,  i.  e, ,  they  are  crpp'er  toleram. 


problem  v  fH  therefore  be  restricted  to  stationary  targets  and  clatter 


. . .  - 


■ 


Hi  f> 


AoectolBji  to  Bqssioi*  {5-2?-  Vitt^m  I5-14J,  A  is  always  diagonal  for  slBasdud  Sasfs 
’Factors.  Ttf 5  mafflstel  ggeperiy  is  fast  a&other  m&Mfejtf  aHto  m  iiw  for*  «h?»  fftpre  *« 
to  essfgf  teffisfer  taUFses  atoiEolsal  fasts  functions  whets  they  az«  ttwingh  a 

tar  syafta*. 


since  A  Is  a  <5s 


n«d  A  A 


seen  nsirices 


I 


(5- 


9  g 

i  eigenvalues  nssoviaic-cl  with  these  eigenvectors  are  j  a  I “  j  a^  |",  '  a^.J 
since 


_  *  T 
[A  A 


XI  3  a  =  0 


or 


O 


O 


;mto:‘orr.:lr-*!cr.  fejction*  (ffl  is  easily  petfarted  fcy  motion  af  the  es?l  - *r.rr.:-ot  a-  to 

“t*c£/'.'cr,  liti)  Is  fiw  sensiffce  to  sHjfet  ttaa^r  fo&cfion  cr.nr.ge5.  arj  (iv  dors  or' 

ferorwe  my  iadicailcB  tot  to  rtangBS  toe  accraed.  tt  stoafcl  ate©  to  Mtftoo  Ort, 

when  mtiMpatli  effects  are  tlitwi»ce4  &  imtfMitsaal  esperleocea  acre  «erf«s 

20 
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5,5  —Continued, 


(0) 


j  U  (a*)  F  (n>) 


da? 


1 

2  JT 


E  (oj)  I  dco 


(5-30) 


where  E  (a?)  =  F  (w)  U  (w)  is  the  t  cho  spectrum. 


Just  <i8  in  the  point  target  case,  maximum  correlator  response  is  pro¬ 
portional  to  echo  energy.  It  is  therefore  relevant  to  derive  a  signal  u  <t>  that  maximizes 
the  reflected  energy  from  a  given  target. 


Writing  (5-30)  in  terras  of  u  (t)  and  changing  the  order  of  integration 


34, 3a 
gives  ' 


/' 


CO  co 


■  ■ 


e  (t)  |  dt  -  /  [  u  (x)  u  (y) 

— Zo 


~  f  |FM!2aJ“'(y-X)da> 


CO  CD 


I  I  ■■ 


(x)  u  (y)  (y  -  x)  dxdy 


(5-31; 


-CC  -C5 


where  R  (r)  is  the  autocorrelation  function  of  the  target  impulse  response  f  (t),  1,  e. , 


— Continued, 


VT) 


f(t)f  (l+r)dt 


(5-32) 


Letting 


gives 


(v)  - 


S3 

/  “ 


(x)  Rff{y-x)dx 


(5-33) 


1  I 


1 


!  e  (t)  !  dt  = 


/  u’iy: 


(y)  2  (y)  dy 


(5-34) 


The  left  hand  side  of  (5-31)  and  (5-34)  is  the  energy  of  the  echo  waveform,  e  (t).  As¬ 
suming  that  both  z  (y)  and  u  (y)  are  normalized  to  contain  unit  energy,  (5-34)  can  be 
maximized  via  the  Schwarz  inequality.  The  right  hand  side  of  (5-34)  is  maximized 
when 


z  (y)  =  K  u  (y) 


Xu(y)  = 


where  >.  is  constant. 


m 

f  U  (s)  Rff  (y  -  X)  d  X 


(5-35) 


(5-36) 


amrmmmmmmmH  mu . . .  ■ r-nTMWEBI 
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5. 5  — Continued. 

If  u  (t)  is  time  limited  to  the  interval  [0.  T  ],  then  (5-36)  becomes  a  homo¬ 
geneous  Fredholm  integral  equation  whose  solutions  are  the  eigenfunctions  of  FL*  (r). 

The  signal  that  maximizes  the  echo  energy  for  a  particular  target  is  the 

solution  of  (5-36)  with  maximum  eigenvalue-  For  this  signal,  the  echo  energy 

is  X  E  . 
max  u 

If  the  transmitted  signal  is  not  theoretically  time  limited  (for  example,  if 
u  (t)  is  a  CW  sipial  with  Gaussian  envelope)  then  the  Fourier  transform  of  (5-38)  is 

X  U  <oj)  =  U  <u?)  [  F  (w)  { 2  (5-37) 

If  U  (w)  is  sufficiently  narrowband,  then 

X  U  (&)  «  j  F  (u?o)  j 2  U  (a?)  (5-38) 

where  is  the  carrier  frequency  cf  u  (cj).  The  narrowband  signal's  carrier  frequency 
should  therefore  correspond  to  a  maximum  on  the  curve  of  baekseattered  power  versus 
frequency.  Such  a  waveform  is  optimum  for  the  maximization  of  echo  energy. 

For  a  point  target, 

R.£  (y  -  x)  =  5  (y  -  x)  (5-38) 


5-21 


E 


il 

mm 
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5,  5 _  —  Continued, 


and  (5-26)  is  identically  satisfied  for  all  admissible  waveforms.  A  planar  target  does 
not  favor  one  signal  over  another  as  far  as  returned  energy  is  concerned.  For  a  non- 
planar  target,  however,  both  echo  energy  and  maximum  echo  power  depend  upon  the 
transmitted  waveform. 


Equation  (5-38)  is  equivalent  to  Equation  (5-29),  where  U  (u)  is  described 
by  one  of  the  eigenvectors  in  (5-28).  The  optimum  eigenvector  is  the  one  with  the 
maximum  eigenvalue,  X  «  Just  as 

!H3L?C 


x 

max 


9 

max  |  F  (w  )  | 
a*.  '  o  1 

r* 


(5-40) 


Ij 


?  ? 


in  (5-38),  Equation  (5-29)  implies  that 


A 


max 


max  (  !  a, „  ! 


i'5-41) 


For  a  sinusoidal  basis. 


3 

m 


9 


(5-42) 


The  above  equivalence  Is  not  surprising.  In  the  absence  of  clatter,  the 

SIR  is  maximized  by  maximizing^  ,  with  constraints  upon  signal  energy,  b  |u  I 

*  it®  2  n  n 

and  upon  the  energy  of  the  filter’s  impulse  response,  b  \  ^  |  .  Since  ,  the  magni¬ 
tude  of  the  filter’s  response  to  the  target  echo,  is  bounded  by  the  energy  of  the  echo, 


f; 


5-22 
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5,  5  — Continued, 


the  SIR  maximization  problem  is  identical  to  maximizing  the  quantity  in  (5-  30}  with 
equivalent  energy  constraints. 

The  equivalence  between  SIR  maximization  and  maximization  of  echo  energy 

shows  that  the  eigenfunction  equation  (5-16)  is  analogous  to  Equation  (5-36),  It  follows 
T  * 

that  the  matrix  A  A  is  analogous  to  (y  **  x).  Indeed,  £  (t)  is  the  target's  impulse 
response,  the  function  that  is  convolved  with  the  signal  u  (t)  in  order  to  form  the  echo. 
Similarly,  the  vector  A  u  describes  the  first  N  terms  of  a  convolution  process,  as  Indi¬ 
cated  by  Equations  (4-19)  -  (4-24). 

One  can  further  exploit  the  equivalence  by  noting  that  the  general  SIR  prob¬ 
lem  defines  A  as  in  (5-13),  Each  element  In  (5-13)  involves  the  difference  between 
target  and  clutter  components.  The  analogous  result  is  to  rewrite  (5-36)  as 

te 

a  u  (y)  =  /  u  (x)  H,  (y  -  x)  d  x  (5-43) 

j  f-C,  I  -  c 

where  R  .  (t)  is  the  autocorrelation  function  corresponding  to  F  (w)  -  C  (m,  the 

I  ■*  O*  I  c 

difference  between  target  and  clutter  transfer  functions, 

5. 6 _ Summary  of  Section  5. 


The  use  of  certain  wideband  signals  allows  for  the  elimination  of  nonlinear 
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of  SIR  becomes  much  simpler.  The  determination  of  an  optimum  signal-filter  pair,  how¬ 
ever,  still  depends  upon  the  basis  that  is  chosen  to  describe  the  relevant  functions.  A 
basis  composed  of  the  functions  (2-S)  seems  to  produce  much  more  practical  solutions 
than  a  basis  of  narrowband  sinusoids. 

The  solution  to  the  SIR  problem  is  similar  to  the  solution  of  a  related  prob¬ 
lem.  !,  e. .  maximising  the  energy  of  the  target  echo.  The  energy  maximization  problem 
becomes  identical  to  the  SIR  problem  if  the  target  transfer  function  F{*>)  is  replaced  by 
P(oj)  -  C  (a?),  the  difference  between  target  and  clutter  transfer  functions.  Both  prob¬ 
lems  reduce  to  the  solution  of  an  eigenfunction  equation.  This  result  suggests  that  the 
functions  (2-8)  could  prove  useful  in  providing  alternative  solutions  to  such  equations. 

In  order  to  maximize  SIR,  the  approaches  in  Sections  2. 3  and  5. 3  both  maxi¬ 
mize  the  difference  between  target  and  clutter  responses,  while  constraining  the  response 
to  noise.  For  square  law  detection,  uniformly  distributed  clutter  has  the  same  effect  as 
noise,  i.e. ,  undesired  responses  occuring  at  random  times.  Clutter  is  thus  treated  as 
though  it  were  noise.  For  linear  detection,  the  difference  between  target  and  clutter 
responses  is  the  same  as  the  response  of  the  linear  filter  to  a  signal  that  is  the  difference 
between  target  and  clutter  echoes.  The  optimum  rece  iver  therefore  acts  as  a  matched 
filter  —  matched  to  the  difference  between  signal  and  expected  clutter  returns. 

The  above  observations  are  easily  applied  to  the  receiver  m  Figure  4-1.  ff 
square  law  detection  is  used,  one  should  first  determine  the  expected  output  power  of  each 
filter  with  noise  and  clutter  echoes  at  the  input,  m  order  to  whiten  the  response  to  inter¬ 
ference,  each  energy  normalized  filter  output  should  then  be  divided  by  its  own  measured 


. . . . . . . . . . . . . . . . mm . . 


ESI?- PEI  15 


5, 6  —Continued. 


response  to  noise  *,nd  clutter.  (If  the  interference  has  constant  power  spectral  density, 

all  energy  normalized  filter  outputs  will  be  divided  *>v  the  same  number. )  This  procedure 

results  in  a  correlation  of  echoes  with  f  /  I  (w) ,  where  I  (w)  is  the  response  to  inter- 

n  A 

ference,  as  in  Eouation  2-35.  For  linear  detection,  one  should  replace  the  gains  f  in 
*  »  n 

Fif  y\.  -i-l  by  f  -  c  ,  where  the  f  coefficients  correspond  to  the  target  transfer 
n  *  ti  n 

function  and  the  coefficients  correspond  to  the  average  clutter  transfer  function. 


ESL-PR1 


G. 


BAT  WAVEFORMS  AND  THE  SIMILARITIES  BETWEEN  AMPLITUDE 
MODULATION  AND  LINEAR  FILTERING  FOR  CHIRPED  PULSES. 


6. 1 _ Bat  Signals  and  Amplitude  Modulation. 


The  experimental  work  of  Griffin,  Friend, 


34, , 


seems  to 


suggest  that  the  tat  Myotis  lucifugus  may  use  amplitude  modulation  differences  to  dis¬ 
tinguish  between  different  target  shapes.  If  the  theory  presented  in  Section  2  is  applied 
to  products  in  the  time  domain  (rather  thar  in  the  frequency  domain)  one  arrives  at  an 
optimum  time  function  (t)  for  Taylor  seties  characterization  of  the  modulating  sig¬ 
nal  a  (t).  This  function,  U7  (t),  is  identic;  1  to  the  function  in  Equation  (2-8),  provided 
the  independent  variable  is  time  ratter  than  frequency. 


it  has  been  found  that  the  signals  U 


(t)  can  be  matched  to  tat  signals  if 


U .  (t)  is  pro  multiplied  by  a  truncating  function  (t),  where 


n=M 


(6-1) 


When  the  signal  (t)  (t)  is  used  for  target  characterization,  the  first  M  coefficients 

of  the  modulation  a  (t)  are  unaffected  by  the  truncation.  The  truncating  function  that 
has  been  used  is 


exp  t  -  0H)  '*  1 


n=l 


(6-2) 


6-1 
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6.1  —Continued. 


Some  examples  of  the  functions  (t)  y,  (t)  are  shown  in  Table  6-1.  The  corresponding 

signals  are  shown  tn  Figures  6-1  and  6-  1.  These  signals  compare  favorably  with  the 

30 

Myotls  erasing  and  pursuit  pulses  shown  in  Figures  6-3  and  5-4. 

G.  2  The  Fourier  Transforms  of  Bat-  Like  Signals. 


The  analysis  in  the  Appendix;  indicates  that  the  inverse  Fourier  transforms 
of  the  functions  U  (_•)  in  (3-S)  have  approximately  the  same  mathematical  form  as  the 
frequency  domain  functions.  This  similarity  implies  that  any  given  signal  (t)  can  be 
approximated  by  the  inverse  Fourier  transform  of  U  («),  where  U  (*)  and  U  {•) 

**  i  w 

belong  to  the  class  of  functions  described  by  Equation  (2-8). 


Figures  6-1  through  6-4  demonstrate  that  some  bat  signals  can  be  matched 

to  the  functions  U,  (ti  v  (t).  The  similaritv  between  U  (t)  and  the  inverse  Fourier 
1  '  1  ’  1 

transform  of  U_  (a?)  would  seem  to  imply  that  the  same  bat  signals  can  be  matched  to 
the  Inverse  Fourier  transform  of  U9  (w)  multiplied  by  an  appropriate  truncating  func¬ 
tion,  y  (cr).  Such  matches  can  indeed  be  made,  provided 

y9  (w)  -  exp  [  -ip  ic) 

In  (6-3),  y  (u:)  acts  as  a  high  pass  function.  The  first  M  Taylor  series  coefficients  cor¬ 
responding  to  F  (w)  y  (oj)  will  be  those  of  F  <w)  alone  if  F  (w)  has  no  more  than  M 
siaimcant  terms,  since 


1  4- 


♦  >: 


nl  0 


I 

1 

I 

I 

1 

i 

I 

I 

1 
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Table  6-1.  Examples  of  the  Signals  (t)  y1  (t) ,  where  U.  (t)  is  Gives  by  EquVion 
(2-8)  sad  yi  (t)  is  an  Appropriate  Truncation  Function.  When  ilr  se 
Waveforms  are  Amplitude  Modulated,  they  Convey  Inform  at  T;.-  about  the 
Taylor  Series  Ccefflcieais  of  the  Modulating  Time  Function 


1. 


U,  (t)  ¥  (t)  =  t 


5e-(xogt) 


COS  L  400  “  {log  t)  i  e 


-{3t)‘ 


UjWy^t) 


30 

-1. 2  5  (log  l)2  cos  C 160  r  flog  I)  3  e~(^ 

t  e 


6-3 


My  oil  p  luoifugna .  M'oaaurod  by  D.  A.  CaMfliMU'i? 


rcBL-  punr> 


0-7 


Pursuit  Pulse  of  Myotis  lucifugus.  Measured  by  D.  A.  Cah  lander  and  J.  J.  G.  Me  Cue  at 
MIT  Lincoln  Laboratory. 
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G.  2  —Continued. 


and  the  first  M  terms  of  a  power  series  e;<pansion  of  F  (w)  y  (a*)  are  identical  with 
those  of  F  («).  even  if  condition  (6-5)  does  not  hold  true. 

The  functions  (a?)  y?  (cu)  whose  inverse  Fourier  transforms  match  the 

Mvotis  pulses  are  shown  in  Table  6-2,  and  the  inverse  Fourier  transforms  are  illus- 

M 

trated  in  Figures  6-5  and  S-6.  The  value  of  P‘  for  the  second  function  m  Table  6-2 
(the  pursuit  pulse)  indicates  that  the  function  in  Figure  6-6  is  capable  of  nearly  un- 
di  snorted  characterization  of  all  the  power  series  coefficients  f  corresponding  to  the 

il 

target  transfer  function,  F  (w). 


6^3 _ Discussion  of  the  Mvotis  Results. 

Figures  6-1  through  6-6  demonstrate  that  two  different  hypotheses  can  be 

used  to  explain  the  Myotts  erasing  and  pursuit  pulses.  One  hypothesis  is  based  upon 

the  observation  that  different  targets  cause  the  echo  of  a  chirped  pulse  to  acquire,  dif— 

35 

ferent  amplitude  modulations.  The  other  hypothesis  is  based  upon  the  idea  that 
targets  act  as  linear  filters. 

These  hypotheses  lead  to  similar  signal  design  problems.  The  problems 
are  similar  because  a  duality  exists  between  (i)  amplitude  modulation,  which  multiplies 
the  transmitted  time  signal  u  (i)  by  a  time  function  a  It),  and  (ii)  linear  filtering  ,  which 
multiplies  the  spectrum  U  {«)  of  the  transmitted  signal  bv  a  frequency  function  F  (tu). 

One  would  expect  the  solutions  of  the  two  problems  to  be  dissimilar,  since 
comparatively  few  functions  are  identical  in  time  and  frequency  (i.  e. ,  only  a  restricted 


6-9 
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Taole  6-2,  Examples  of  the  Transforms  U  (tu)  y  {&),  wise  re  U  (a;-)  is  Given  by 

Equation  (2-8)  and  y2  (co)  is  an  Appropriate  Tnin eating  Fui:ction,  When 
the  Inverse  Fourier  Transforms  of  U_  {&>)  y„  (w)  are  Passed  Through  a 

M  ^ 

Linear  Filter,  they  Convey  Information  about  the  Taylor  Series  Coef¬ 
ficients  of  the  Filter’s  Transfer  Function, 


2  -20 

T7  .  ,  ,  .  -0. 7  -0.  S  Ci og  co)  -j  241.  2  r  lo^  u  -{.  9a0 

U2  <w)  y2  {<u)  =0;  e  ■  e  'o' 


2  -22 

U,  M  v  (M)  =  „-»• 1  e-°- 6  (i0g  ">  <T*  80  T  l0*  "  e-»- 8  "> 


ft  Ml 


The  Inverse  Fourier  Transform  of  the  First  Function  in  Table  G-2. 
Like  the  Porpoise  Waveforms  in  Figure  2 -(>,  This  Signal  is  Designed 
for  Characterization  of  a  Linear  Filter  Via  Spectral  Power  Series 
Coefficients.  Instantaneous  Period  Incroases  Linearly  with  Time. 
Compare  with  Figures  G-.l  and  6-3. 
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class  of  functions  are  exact  eigemalues  of  the  Fourier  transform).  Chirped  signals, 
however,  are  often  approximate  eigenvalues  of  the  Fourier  transform.  Linear  chirp 
is  the  best  known  case,  but  the  stationary  phase  approach  (see  the  Appendix)  implies 
that  some  other  chirped  waveforms  have  the  same  property. 

Given  that  two  mathematical  hypotheses  result  in  approximately  the  same 
correspondence  with  reality,  one  must  seek  further  data  to  determine  which  hypothesis 
is  correct,  or  whether  there  exists  an  equivalence  of  the  hypotheses  that  has  been  over¬ 
looked. 


The  hypothesis  that  bat  signals  are  designed  to  analyze  amplitude  modu¬ 
lation  is  re-enforced  by  some  neurophysic  logical  work.  In  studies  of  bats’  auditory 
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neural  processing,  N.  Suga  and  A.  GrinneH**  have  found  many  neurons  which  are 
tuned  to  particular  frequency  intervals.  Since  neuronal  excitation  is  a  function  of  amplitude 
as  well  as  frequency,  each  tuned  neuron  can  act  as  an  AM  demodulator  over  a  particular 
frequency  range.  Tuned  neurons  therefore  provide  a  means  for  determining  the  echo 
amplitude  at  any  given  frequency. 

The  hypothesis  that  bat  signals  are  designed  to  convey  information  about 
the  target  qua  linear  filter  is  strengthened  by  the  observations  in  Sections  2  through 
i.  The  psychoacoustic  excitation  functions  can  theoretically  be  inv  oked  to  explain  not 
only  porpoise  echolocation  clicks,  out  those  of  any  ether  mammal, 

A  supporter  of  the  amplitude  modulation  hypothesis  would  have  to  concede 
that  there  is  considerable  evidence  to  suppoit  the  critical  bandwidth  theory.  Neverthe¬ 
less,  he  could  point  out  that  N.  Suga  has  found  no  evidence  of  a  dispersive  delay  pulse 
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compression  mechanism  in  his  years  of  studying  the  FM  bats.  One  could  therefore 

claim  that  the  auditorv  filters  do  not  implement  the  log  phase  necessary  for  true  pulse 

40 

compression  and  efficient  linear  filter  characterization. 

It  Is  possible,  however,  that  FM  bats  may  use  a  pulse  compression  mecha¬ 
nism  that  does  not  inv  olve  a  dispersive  delay.  Fez  example,  incoming  signals  may  be 

% 

repeatedly  multiplied  bv  stored  versions  oi  the  transmitted  waveform  ,  and  the  product 
integrated.  The  multiplication  operation  is  particularly  easy  to  perform  with  neural 
elements,  sine~  a’l  signals  are  cooed  In  terms  of  a  positive  quantity,  the  number  of 
excitation  pulses  per  second.  Positive  voltages  can  be  interpreted  as  increased  dis¬ 
charge  rates  (excitations)  and  negative  vol  ages  can  he  interpreted  as  decreased  rates 
(inhibitions),  relative  to  the  spontaneous  discharge  rate  of  the  unexcited  neuron.  In 
electronics,  the  multiplication  of  two  j.osit  ve  voltages  (single  quadrant  multiplication) 
is  accomplished  with  a  voltage  variable  amplifier.  Any  neuron  with  both  excitatory  and 
inhibiiorv  synaptic  connections  to  its  dendi  -tes  provides  the  analog  of  such  an  amplifier. 
The  integration  process  can  be  implement.  J  by  low  pass  filtering,  which  again  is  natu- 
callv  realized  by  any  neuron;  all  neurons  hive  a  recovery  time  that  limits  the  maximum 
discharge  rate,  so  that  they  act  as  low  pass  filters. 

A  supporter  of  the  hypothesis  that  targets  are  characterized  as  linear 

filters  can  therefore  still  insist  that  the  apircpriate  filtering  process  (Figure  4-1)  is 

implemented  by  Mvotis.  He  can  point  out  that  only  ioo  neurons  are  required  for  the 

*  * ** 

whole  correlation  process.  Furthermore,  he  can  remark  that  the  reference  signal 


*  For  a  20  KHz  bandwidth,  this  method  would  require  a  new  multiplication  to  commence 
every  50  microseconds, 

**  Any  one  correlator  can  be  reused  after  t  has  performed  multiplication  and  integration 
over  the  expected  signal  duration,  l  or  a  ;  msec,  20  kHz  wide  signal,  the  receiver  could 
thus  be  implemented  with  100  multiplicatio  i- integration  circuits. 

e-i4 
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is  based  upon  an  efferent  input  from  a  higher  part  of  the  neural  processor,  i,  e. ,  a  part 
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of  the  brain  that  could  be  incapacitated  by  anesthesia.  These  observations  could  explain 
the  lack  of  success  that  has  been  encountered  in  searching  for  a  pulse  compression  mecha¬ 
nism  via  neural  probes  on  anesthetized  animals. 


The  Equivalence  of  the  Amplitude  Modulation  and  Linear  Filtering 
Hypotheses. 


Since  additional  evidence  can  be  found  to  favor  either  hypothesis,  let  us 
reconsider  the  possibility  that  the  hypotheses  arc  equivalent.  As  it  turns  out,  there 
is  considerable  evidence  favoring  this  point  of  view: 

First,  there  is  the  fact  that  if  the  chirp  rate  is  sufficiently  slow,  the  FM 
signal  determines  the  response  of  a  linear  filter  to  each  frequency,  and  this  response 
is  manifested  as  an  amplitude  modulation  of  the  chirped  waveform. 


Second,  the  amplitude  of  a  chirped  signal  in  time  is  easily  related  to  its 
amplitude  in  frequency.  If  the  stationary  phase  approximation  is  valid  (see  the  Appendix.) 
It  follows  that  a  giver,  amplitude  modulation  in  time  results  in  a  corresponding  amplitude 
modulation  in  frequency  (linear  filtering),  i.  e. ,  there  is  an  equivalent  linear  filter  for 
anv  amplitude  modulation  imparted  to  the  chirped  signal. 

Third,  the  matched  filtering  shown  in  Figure  4-1  results  in  a  process  that 
is  similar  to  sampling  the  magnitude  of  an  echo  over  different  segments  of  the  returned 
waveform.  The  signals  u*  (knt)  are  compressed  when  n  is  positive,  so  that  echoes  are 


— Continued, 
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correlated  with  functions  whose  maxima  occur  closer  and  closer  to  the  beginning  of  the 
pulse  as  n  increases.  The  same  signals  are  stretched  when  n  Is  negative,  so  that  echoes 
are  correlated  with  functions  whose  maxima  occur  further  and  further  toward  the  end  of 

A 

the  pulse  as  n  decreases.  This  effect  wiU  be  demonstrated  in  Section  7.  The  effect 
may  also  inmeate  a  scheme  by  which  the  tuned  neurons  found  by  Suga  and  GrinneU  can- 
act  as  the  equivalent  of  a  tank  of  matched  filters  or  correlators. 


Fourth,  the  two  derivations  of  Mvotis  signals  hinge  upon  the  use  of  power 
series  to  describe  the  amplitude  modulation  function  a  ft)  or  the  linear  filter  function 
F  (*c).  It  has  already  teen  shown  (Section  4.  4)  that  linear  filtering  results  in  a  convo¬ 
lution  of  filter  coefficients  with  the  spectral  coefficients  of  the  processed  signal.  In 
fact,  amplitude  modulation  can  also  be  written  as  a  convolution  operation  involving  the 
spectral  coefficients  of  the  processed  signal.  This  effect  is  discussed  below. 


6.  5  A  Mathematical  Relation  between  Amplitude  Modulation  and  linear 

Filtering. 


interval  (0  , 


Consider  the  product  of  two  functions  a 

_  -  - 

T)  .  a  (t)  can  be  represented  by  an  N 


(t)  and  u  (t).  On  a  finite  time 

ft 

order  polynomial. 


a(t)  =  Ao 


(6-8) 


" - - — -  - -  n  *  ^  n 

*  The  signals  (3-8)  have  no  range- doppler  coupling,  so  that  u"(k"t)  * — >  U®  (u?/  k  )  /  k  ' 
correlates  with  u  (t)  for  zero  delay.  Zero  delay  implies  that  u*(knt)  and  n  (t)  have  the 
same  time  origin  (t  =0>  for  maximum  correlator  response.  Plots  oi  u*  (knl)  for  different 
values  of  n  reveal  maxima  whose  locations  depend  upon  n. 

’  *  As  in  (6- 1),  we  represent  power  series  coefficients  of  time  domain  ftmctions  by  capi¬ 
tal  letters  and  coefficients  of  frequency  domain  functions  by  lower  case  letters. 


where  the  aouoie  snw  indicates  a  Fourier  Iran 

power  series  coefficients  c  of  C  to) . 

a  '  7 

According  to  (3-3),  one  way  to  fled 

I™  c  (t)  d  t  =  (ml 
o 


Substituting  (6- a)  Into  C6-S)  and  (6-9)  Into  (6- Hi 
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Recall  that  the  coefficients  e_  correspond  to  the  spectrum  C  (10),  where 
C  (w)  is  the  Fourier  transform  of  a  (t)  u  (t).  The  coefficients  c.  can  be  estimated  with 
the  processor  in  Figure  4-4. 

If  C  (u>)  results  from  the  linear  filtering  process 

C  (oj)  =  F(w)  U<w)  ,  (G-  23) 

then  the  power  series  spectral  coefficients  of  F  (w)  can  be  derived  by  computing  T  1  £ . 
where  is  the  NxN  matrix  in  (4-33)  . 

If  C  (cu)  results  from  amplitude  modulation  of  the  signal  u  (t),  i.  e. ,  if 

C  (cv)  < — ♦  c  (t)  =  a  (t)  u  (t)  < — ♦  A  (u>)  *  U  (to)  .  (G-  24) 

then  the  power  series  temporal  coefficients  of  a  (t)  can  be  derived  by  computing  M  1  c  , 

where  M  is  the  NxN  matrix  in  Equation  (6-20). 

—  u 

The  r-bove  two  processes  are  shown  in  Hgure  6-7.  The  first  process 
gives  the  viansfer  function  coefficients  corresponding  to  a  filtered  version  of  the  trans¬ 
mitted  signal  u  (t).  The  second  process  gives  the  equivalent  amplitude  modulation 
coefficients  for  the  same  received  signal. 


G-23 
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6.  ?  _ Estimation  of  Amplitude  .  atton  Parameters  or>  the  Basis  of 

More  Than  One  Ihiise, 


When  a  pulse  train  is  transmitted,  different  received  pulses  will  generally 
have  different  amplitude  modulating  functions.  These  functions  can  be  related,  however, 
since  they  are  presumably  determined  by  the  same  continuous  amplitude  modulation,  a  (t). 
If  the  interpuJse  interval  is  Tq  and  the  modulation  of  the  first  pulse  is  given  by  (6-8),  then 
the  modulation  of  the  second  pulse  is  given  by 


a  { 


T  ) 
o 


N 

2 

n=0 


f  a*n)  (T  )  /  n  !  J  (  t  -  T  Y 
o  o 


(6-25) 


yn) 

where  a'  (T  )  is  determined  bv  differentiating  both  sides  of  (6-8)  and  letting  t  -  T  ; 

o  o 


•w  <V 


y 


(m  1  A  )  T  m  "  n  /  (m  -  n)  ! 
m  o 


(6-  28' 


m-n 


Substituting  (6-26)  into  (6-25), 


N  (  t  -  T  )”  N 
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ml  .  _  m - n 

- 1 

(m  -  n)  I  m  o 


(6-28) 


In  matrix  notation, 


!  1  T 


o  o 


N(N-l)  T 


A  1 

N1  i 


O 


Av  1 

N  l 


Defining 


a  ( t  -  2T  )  =  Y  A  (  t  -  2T  )f 

o'  n2  c 


(6-29) 


Y  C  a^n>  (2T  )  /  n  •  3  (  t  -  2T  )n  , 


(6-30) 
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we 


can  determine  a^  (  2T  )  by  differentiating  both  sides  of  (6-  27)  and  letting  t  -  3T  ; 


N 


a'"'  (2Tq)  =  ^ 

m=n 


m!  A  T  m  n  /  (m  -  n)  I 
m  1  o 


(6-31) 


Substituting  (6-31)  into  (6-301, 


a(l-2To) 


|  A„2 

n-0 


N  „  K 

2  ‘n-  2 

n=0  m=n 


mi  _  m  -  n  .  _  n 

- - —  T  A  j  (  t  -2T  ) 

(m  -  n)  I  o  ml  o 


(6-32) 

so  that 


Ao2 

1 

'P 

» 

o 

T  2  . 

0 

• 

* 

_  N 

T 

0 

A 

o  1 

A12 

0 

1 

2T 

0 

NT  N'1 
o 

A 

i  i 

1  ! 

' 

* 

i : 

A1 1 

A 

n 

n 

1 

N-2 

K(N-l)  T 

0 

A 

A2  2 

V 

2  1 

2  J 

rt  . 

21 

4 

• 

• 

4 

• 

• 

4 

• 

- 

AN2 

o 

0 

0 

• 

• 

1 

AN1 

— 

_ 

6- 
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(6-36)  , 


r  ■an  I-1 
u  P  j 


(6-37) 


A  ! 

Nn  |  (6-38) 

is  the  vector  of  *<me  series  coefficients  measured  from  the  n—  received  pulse. 


The  coefficients  measured  for  alMhe  received  pulses  can,  mutatis  mutandis, 

be  used  to  estimate  the  time  series  coefficients  in  (6-8)  .  Each  pulse  is  processed 

as  in  Figure  6-7  (the  second  process).  The  resulting  coefficients  can  then  be  used  to 

estimate  A  ,  A  .....  A^.  in  (5-8)  by  means  of  the  transformation  (6-36). 
o  l 


. . . . . * 
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6.  S _ Discussion  and  Summar>/  of  Section  6. 


The  theo^’  of  target  description  that  has  been  derived  to  explain  cetacean 

echolocation  has  been  applied  io  bats.  The  theoretical  signals  can  be  matched  to  Myotis 

* 

pulses  as  well  as  to  those  of  Tursiops,  To  match  the  bat  pulses,  the  theoretical  wave¬ 
forms  must  be  high  pass  filtered  in  a  special  way.  The  high  pass  filter  must  not  dis¬ 
turb  the  ability  to  determine  at  least  twenty  Taylor  coefficients  of  a  target's  transfer 
function. 


The  above  results  constitute  the  first  mathematical  evidence  that  a  unified 
theory  of  animal  echolocation  exists.  Tlie  theory  also  provides  a  measure  of  the  maxi¬ 
mum  number  of  spectral  coefficients  that  can  convenient]}  be  estimated  b\  the  animals. 

Because  of  similarities  between  the  theoretical  waveforms  and  their 
Fourier  transforms,  one  can  interpret  the  tat  signals  (and  possibh  the  porpoise  signals 
as  well)  in  terms  of  lime  domain  charaeieri/.ation  of  an  unknown  amplit  ide  modulation. 
This  alternate  interpretation  is  interesting  from  several  viewpoints: 

(i)  It  {mplies  that,  for  certain  classes  of  signals,  one  can 
find  a  time  invariant,  linear  filtering  operation  that  is 
equivalent  to  any  amplitude  modulation.  This  equivalence 
implies  (Section  10)  the  existence  of  a  time  invariant  filter 
which  has  the  same  effect  as  any  time  varying  linear  trans¬ 
formation  of  a  given  pulse. 


*  ihe  theoretical  signals  have  also  been  matched  to  echolocation  pulses  of  Ep.esicus 
fuscus  and  Lasiurus  borealis.  These  matches  will  be  discussed  in  a  future  report. 


6-30 
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(ii)  The  AM  interpretation  provides  insights  Into  the  possible 
functions  of  various  neurological  processes  which  have 
been  discovered  in  the  auditory  pathways  of  bats  and  other 
animals.  Specifically,  one  can  interpret  the  functions  of 
neurons  which  respond  to  specific  frequencies  or  chirp 
rates. 

(iii)  The  AM  interpretation  illustrates  another  virtue  of  power 
series  spectral  analysis,  viz. ,  that  both  multiplication 
and  convolution  In  the  time  domain  can  be  represented  by 
the  same  operation  (convolution  of  appropriate  coefficients) 
in  the  frequency  domain. 

in  a  discussion  of  neural  processing  (as  seen  from  the  viewpoints  o? 
amplitude  demodulation  vs.  linear  filter  characterization)  it  was  pointed  out  that  a 
single  matched  filter  could  be  implemented  with  100  neurons,  orovided  suitable  inputs 
were  available.  It  would  seem  that  the  proper  input  signals  can  be  made  available 
via  the  parallel  connection  of  neurons  to  carry  high  frequency  information.  This 
concept  is  discussed  further  in  Section  8. 

It  is  possible  that  correlation  using  a  single  nerve  cell  could  be  performed 

at  a  very  peripheral  level,  perhaps  at  the  location  where  pressure  waves  are  first 

converted  into  chemical-electrical  signals.  Such  a  hypothesis  would  require  (i)  that 

41 

optimum  filtering  be  detectable  at  the  peripheral  level  if  anesthetic  Is  not  used  , 


6-31 
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6.  8 


(ii)  <hat  efferent  ner  e  fibers  extend  all  the  way  to  the  vicinity  of  the  hair  cells''4  ,  and 
(HI)  that  the  correlation  processing  of  a  wideband  signal  be  divided  into  a  sequence  of 
narrow  nan  a  correlations  over  restricted  frequency  ranges.  * 3 
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AN  EXPERIMENT  TO  TEST  THE  THEORY. 


1 _ The  Physical  Implementation. 

In  order  to  gam  further  insight  into  tlie  iheorj  and  how  it  works  in  practice, 

■* 

a  sonar  experiment  is  being  carried  out  with  high  frequency  sound.  Although  the  ex¬ 
periment  was  funded  by  ESL  as  a  x eduction  to  practice  for  a  patent  application,  it  is 
reported  here  in  order  to  augment  the  reader's  comprehension  of  the  theory. 


The  ultrasonic  transducers  that  were  used  for  the  experiment  were  built 
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to  the  specifications  given  by  J..J.  C-.  Me  Cue  ami  A.  Bertolini  .  Although  the  trans¬ 


ducers  were  designed  as  receh  ers,  one  was  used  as  a  transmitter.  The  two  trans¬ 
ducers  are  shown  in  a  bistatic  configuration  in  Figure  7- 1.  The  target  is  the  plane  in 
the  foreground;  it  Is  perpendicular  to  the  line  that  bisects  the  angle  between  the  trans¬ 
ducers. 


fhe  signal  that  is  being  used  has  a  Fourier  transform  as  in  Equation  (3-8). 
Specifically,  the  signal  has  Fourier  transform 


U  {, 


.  85  -{log  u?)  /  2iog  K  -  j  16  «  log  «:  /  log  k 


where  k  *  2.  06.  This  value  of  k  is  slightly  higher  than  the  values  usually  employed 
by  bats  and  dolphins.  The  reason  for  such  a  large  value  was  so  obtain  better  isolation 
between  coefficients,  since  we  planned  to  measure  only  six  different  target  coefficients. 


P. 


The  experiment 
Keleshian:  data 


was  designed  by  R.  Alles;  the  electronic  equipment  was  built  by 
is  being  gathered  and  processed  by  W.  Reese. 


KSI,  I  Mil  \:> 
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The  signal  was  chosen  to  be  a  compromise  between  oat  jmd  porpoise  echo- 
location  waveforms.  Like  dolphin  signals,  no  truncation  function  was  applied  in  tue 
frequency  domain.  like  bat  signals,  the  phase  factor  was  made  large  enough  to  pro- 
dace  a  significant  number  of  zero  crossings,  i.e. ,  ft  recognizable  chirp.  The  ptess 
factor,  although  large  enough  to  produce  a  chirp,  is  still  not  as  large  that  employed  by 
Mv&tis-  hut  much  larger  than  that  on: cloyed  by  Torsions. 

The  transmitted  signal  Is  shown  in  Figure  7-.'.  The  waveform  Is  generated 
in  toe  laboratory  by  using  a  read- car/  memory  (ROM)  dfc/toe.  By  varying  the  clock 
rate,  L  e. ,  the  rate  at  which  the  sampled  signal  is  read  out  rf  the  ROM,  the  waveform 
can  be  stretched  or  compressed.  Notice  that  the  signal  in  Fi  7-2  is  truncated  at 
its  low  frequency  end.  This  truncation  was  decided  uccr.  because  the  low  frequency  etui 

has  comparatively  little  energy  and  is  greatly  mer  sampled,  The  samples  could  bet- 

* 

ter  be  utilized  at  the  high  frequency  end  of  the  puls:-. 


The  correlation  operation  is  performed  with  analog  differentiators  aijd  inte¬ 
grators  and  an  .analog  multiplier.  The  reference  signal  is  Stained  by  generating  a 
second  waveform  with  the  ROM,  The  second  waveform  Is  generated  at  the  right  time 
so  that  it  coincides  with  the  target  echo.  The  reference  signal  is  then  differentiated 
once,  twice,  or  three  times,  integrated  once  or  twice,  or  left  alone.  The  resulting 
waveforms  are,  theoretically,  time  scaled  version:  of  the  transmitted  signal.  The 
extent  to  which  differentiation  or  integration  reproduces  a  scaled  version  of  the  signal 
Is  determined  by  the  similarity  between  the  analog  ©Derations  and  %he!r  mathematical 


•  For  efficient  sampling  of  a  chirped  signal,  one  should  really  use  a  variable  clock 
so  that  many  samples  are  not  wasted  In  oversampling  .he  low  frequency  poriloc  of  tl 
signal. 


1 1  »•,  nyi*  i,<  rtwiiiww  ipii'ir  -Miwwti  'fi  *  ^  r«r  ■  “•hf  11  w  i"i'  Wi'HW 


milled  Signal  for  the  Ultrasonic  Experiment.  The  Signal  is 
-July  Memory)  Device  that  was  Programmed  with  a  Comput 


i 
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7. 1  — Continued. 


eounteroa.  ts.  and  by  the  similarity  of  the  signal  to  its  mathematical  counterpart.  Fig¬ 
ure  7-3  shows  that  two  analog  diffei eniiations  and  one  analog  integration  are  reasonably 
well  behaved.  The  second  integration  is  imperfect  because  of  low  frequency  distortion 
(integrators  are  low  pass  filters)  and  because  the  pulse  is  truncated  for  efficient  ROM 
utilization. 

A  scale!  reference  signal  is  correlated  with  the  target  echo  fcv  multiplication 
and  integration.  The  maximum  output  of  the  correlator  is  determined  for  each  reference 
signal.  The  maximum  outputs  are  related  to  the  target  power  series  coefficients  by  the 
analysis  in  Section  4. 

7.  2  Some  Qualitative  Results. 

< 

I 

l 

The  most  important  results  o  this  experiment  are  the  qualitative  observa¬ 
tions  a.vl  insights  that  are  obtained  before  one  even  begins  to  collect  numerical  data. 

Echo  differences  for  different  objects  are  easily  observed  on  the  oscilloscope 
screen.  These  differences  seem  to  be  manifested  as  different  amplitude  modulations 

3  5 

imparted  to  the  returning  pulse.  This  observation  has  already  been  reported  by  Griffin  . 

Figure  7-3  illustrates  that  the  scaled  signals  have  their  maxima  at  different 
times,  relative  to  the  returning  echo.  When  the  scaled  reference  signals  are  correlated 


**  A  better  (but  more  e.'qpensive)  way  to  produce  the  reference  signals  is  U<  use  a  different 
clock  rate  for  the  second  ROM  pulse. 
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7. 2  --Continued. 

with  an  echo,  the  correlator  output  is  verv  dependent  upon  the  amplitude  of  the  echo  at 
the  peak  of  the  reference  signal.  Matched  filtering  for  power  series  target  charaeteri- 
zation  in  thus  equivalent  to  amplitude  demodulation  of  various  portions  of  the  reflected 
pulse.  This  equivalence  has  been  explored  in  Section  6. 

Further  insight  is  obtained  when  one  realizes  that  the  reference  signals 

which  "look"  at  the  leading  edge  of  the  echo  pulse  are  the  differentiated  or  shortened 

versions  of  the  transmitted  signal.  Differentiated  reference  signals  are  therefore  used 

to  characterize  the  leading  edge  of  an  echo,  which  is  related  to  initial  curvature  a^d 

* 

"hardness”  of  the  material. 

Conversely,  integrated  (stretched)  reference  signals  characterize  the  trail¬ 
ing  part  of  the  echo,  and  are  thus  related  to  the  range  extent  of  the  target.  Experimental 
ly,  targets  that  are  extended  in  range  cause  the  echo  to  be  stretched.  This  stretching 
can  be  characterized  by  an  integration  process  when  the  appropriate  signals  are  used. 


*  T.  R.  Bullock  and  S.  Ridgway  have  found  a  dichotomy  in  Tur slops’  neural  processing 
of  acoustic  data.  Signals  with  fast  rise  times  are  processed  in  a  different  part  of  the 
brain  than  signals  with  slow  rise  times.  Since  transmitted  echoloeation  clicks  generally 
have  fast  rise  times,  Bullock  and  Ridgwav  chose  to  classify  the  fast  rise  time  analyzer 
as  the  sonar  processor.  Slow  rise  time  echoes,  however,  are  still  possible  with  fast 
rise  time  signals.  Such  echoes  would  occur  for  targets  with  gradual  changes  in  acoustic 
impedance  or  cross  section,  such  as  a  large  tilted  plane  or  a  soft,  waterlogged  object, 
and  for  all  targets  that  are  far  away,  since  the  channel  acts  as  a  low  pass  filter.  It  is 
therefore  tempting  to  reclassifj  the  fast  rise  time  analvzer  as  the  "sonar  processor 
for  close,  hard  targets.''  Reflections  from  such  targets  ma>  require  quick,  decisive 
action,  whereas  ether  targets  can  be  contemplated  at  leisure.  This  reaction  time  con¬ 
sideration  may  explain  the  dichotomy  found  by  Bullock  and  Ridgway.  An  interesting 
behavioral  experiment  w-ould  be  the  measurement  of  reaction  time  as  a  function  of  signal 
shape. 
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—  Continued. 


The  author  has  often  beer,  asked  about  the  correlation  of  spectral  power 

series  coefficients  with  target  shape.  The  above  obse n  ations  supply  a  partial  answer. 

Small,  hard  targets  with  limited  range  extent  will  have  comparatively  large  coefficients 

for  n  >  0  ,  where  F  «'uj)  -  V  f  .  Targets  with  large  range  extent  or  slow  changes 
'  -W  2  n 

in  acoustic  impedance  will  have  comparatively  large  coefficients  for  n  <  0  . 


?.  3 


Some  Quantitative  Results. 


A  schematic  diagram  of  the  experiment  is  shown  in  Figure  7-1.  The  second 
ROW  pulse  is  correlated  with  the  echo,  after  being  passed  through  a  differentiator  or 
integrator.  The  maximum  correlator  response  magnitude  determines  a  coefficient 
that  helps  to  characterize  the  target  transfer  function.  The  filter  outputs  should  be 
processed  as  in  Section  4,  Figure  4-3.  For  our  purposes,  however,  it  suffices  to 
show  that  the  correlator  outputs  are  indeed  different  for  different  targets. 


All  correlator  responses  should  be  normalized.  Each  response  should  be 
divided  by  the  output  of  the  autocorrelator,  i.  e, ,  the  correlator  response  when  the  ref¬ 
erence  signal  is  the  same  as  the  transmitted  signal. 


Three  targets  have  bee 


used:  A  steel  cylinder  with  a  3/8  inch  diameter. 


a  metal  plate  oriented  perpendicular  to  the  angle 
penri’cular  plane"),  and  the  same  plane  tilted  th.tr 
(he  transmitter. 


bisecting  the  two  transducers  (the  " par¬ 
ty  degrees,  with  the  closest  edge  toward 
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7.3  — Continued. 


The  data  is  presented  as  reflected  signals  and  numbers  in  Figure  7-5  and 
Table  7-1,  Figure  7-5  shows  the  echoes  from  the  three  diffeicnl  targets.  Table  7-1 
Illustrates  *hat  one  can  indeed  obtain  different  normalized  correlator  responses  for  dif¬ 
ferent  target  shapes. 


Future  Experimentation. 


Some  of  the  data  gathered  in  future  experiments  will  be  digitized 
cessed  by  computer.  Computer  filtering  will  eliminate  the  uncertainties  caus 
imperfections  of  the  analog  equipment. 


and  pro- 
ed  by 


Experiments  are  currently  under  way  to  gainer  data  when  receiver  and 
transmitter  are  colocated,  rather  than  separated  as  in  Figure  7-1.  W c  are  alsc  dis¬ 
criminating  lx? tween  cylinders  of  different  size,  in  an  attempt  to  simulate  some  of  the 
behavioral  experiments  that  have  been  performed  with  cetaceans. 

It  is  hoped  that  the  experiment  will  soon  be  set  up  to  discriminate  targets 
in  water,  using  small,  off-the-shelf  wideband  transducers.  It  would  also  be  interesting 
to  detect  small  disturbances  on  the  air- water  interface,  as  Noctilio  leporinus  (the  fishing 
cat)  apparently  does.  **  *  One  could  set  up  the  capacitive  (air)  transducers  to  look  at  the 
surface  of  an  aquarium  containing  fish,  and  electronically  trigger  a  camera  to  record 
the  fishes'  positions  when  an  appropriate  set  of  filter  responses  is  observed. 
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ri^ure  i-3. 


Graphs  of  the  Echoes  from  Three  Different  Reflectors:  (as  the  3/8" 
Diameter  Cylinder,  (b)  the  Perpendicular  Plane,  (c)  the  Tilled  Plane. 
The  Coefficients  In  Table  7-1  are  Determined  by  Correlating  these 
Echoes  with  the  Sipials  in  Figure  7-3. 
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TW 0  MODELS  OF  THE  AUDITORY  SIGNAL  PROCESSOR. 


Experiments!  Results. 


t 


Measurements  of  audi-ory  threshold  as  a  function  of  tea  duration  have 
been  made  for  dolphins’ 5  and  humane  ’°.  For  all  audible  frequencies,  these  me-sure- 
mcai-s  indicate  as  increase  in  sensitrv  for  pu-sc  durations  up  is  approximately  0.  2 

0.  2  seconds,  tfce  threshold  stavs  aouroximatelv  c%mstant. 


'or  laws  banger  the 


** 


-oa su re ments  of  just  noticeable  frequency  difference  h;\  e  teen  mode  foi 
•17 

human  lirteaers.  "  j  ease  meas-irenKscts  indicate  sensitivity  to  a  3  Hz  frequency  change 
for  tones  with  frequency,  less  than  2  fCHa.  Above  2  KHz*  the  minimum  dsiet 
change  berorees  larger:  10  Hz  at  j-  . »  22  Hz  at  $  KHz. 


able  freqaars 


11 


Masking  of  pure  tones  by  narrowband  noise  has  led  to  a  model  of  the  ear 
as  a  spectrum  analyzer.  This  motel  uses  filters  whose  hancteidths  are  approximately 
proportional  to  their  center  frequencies.  These  so-called  critical  barswidlhs  are  sp- 

,  94 

proxunateh  twenty  times  tne  minim —  ‘  ’  *  — 

above. 


pert-edible  frequency  difference  discussed 


£ 


Physiological  experiments  performed  by  von  Befcesy  and  reported  in 

4h 

Fietcher’s  book '  show  a  relation  between  number  of  nerve  endings  attached  to  the  basilar 
membrane  at  a  given  position  and  the  resonant  frequency  at  that  position  (Figures  S-l  ai^ 
2).  The  curves  imply  that  a  narrowband  masking  signal  offset  by  £  f  from  a  pure  tone 
wiU  succeed  in  masking  the  tone  if  At  is  sufficiently  small.  At  the  higher  frequencies, 

&i  cars  be  increased  and  masking  will  still  occur.  This  observation  is  intimately  related 
to  the  concept  of  critical  bandwidth,  as  indicated  by  the  masking  data  curves  in  Figure  S-2 
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Fletcher’s  Model  of  the  Auditory  Processor 
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R.5.  Gases  has  tainted  cut  that  the  a  be*  e  measaremenis  can  aO  he  ex* 
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--Continued. 


The  combination  of  N  (fQ)  fillers  explains  the  just  noticeable  frequency  dif¬ 
ference  data.  To  translate  aCWtone  from  one  filter  combination  to  the  next,  the  frequency 

must  be  shifted  by  r  N  (f  )  y  5  Hz  1  /  2,  i.  e. ,  half  the  total  bandwidth.  For  f  -s  1  KHz  , 

o  o 

'  N  (f  )  x  5  ]  /  2  =  2.  5  Hz  ;  at  f  =  4  KHz  the  required  translation  is  10  Hz  ;  at  f  = 

0  0  o 

8  KHz  ,  the  required  translation  is  20  Hz.  These  numbers  correlate  closely  with  the 
just  noliceabie  frequency  differences. 


The  tone  duration  data  is  also  explained  by  the  new  model,  since  signals 
are  convolved  with  a  filter  impulse  response  fhat  is  0.  2  seconds  long.  The  maximum 
response  to  an  input  pulse  of  constant  frequency  will  thus  increase  linearly  with  pulse 
duration  until  the  pulse  H  .comes  longer  than  0.  2  seconds.  At  this  point,  the  maximum 
response  remains  constant,  even  when  the  input  pulse  duration  is  increased. 


The  new  model  is  also  consistent  with  the  capacity  of  neurons  to  transmit 

53 

a  signal.  Tue  refractory  period  of  a  tvp’cal  neuron  lasts  approximately  one  millisecond. 

Because  of  this  refractory  period,  a  tone  with  more  than  one  maximum  per  millisecond 

{  1  KHz)  can  only  be  transmitted  if  the  burden  is  shared  by  other  neurons.  Assume  that 

these  neurons  are  connected  to  adjacent  5  Hz  filters.  Spontaneous  discharges  of  different 

neurons  will  cause  them  to  be  excited  at  different  times,  since  the  excitation  of  one  neuron 

may  occur  during  the  refractory  period  of  another.  If  e.eh  neuron  triggers  on  a  different 

signal  peak,  then  a  combination  of  four  neurons  can  handle  a  4  KHz  signal,  eight  neurons 

can  handle  an  8  KHz  signal,  etc.  In  general,  a  combination  of  N  (f  )  neurons  can  handle 

w  o 

a  signal  with  frequency  f  ,  where  N  (f  )  was  defined  above.  If  the  N  (f  )  neurons  are 

o  or  o 

connected  to  adjacent  5  Hz  filters,  the  just  noticeable  frequency  difference  will  be 

r  N  {f  ?  y  5  Hz  1  /  2. 
o 


w  • 


. . .  , 
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8.  3 _  — Continued. 

C.  S.  Johnson  has  pointed  out  that  the  new  model  can  also  explain  the 
rapid  cutoff  of  marine  mammals'  auditory  thresholds  at  the  high  frequency  end.  The 
model  can  also  explain  changes  that  occur  in  the  auditory  spectrum  analyzer  for  w > de¬ 
band  signals.  D.  M.  Green  has  found  that  "critical  bands  are  not  fixed  in  width  but  are 
adjustable  so  as  to  match  the  particular  detection  situation.  ”'J° 

The  most  important  aspect  of  the  new  model  is  that  it  allows  for  faithful 
transmission  o'  '  complete  signal  without  the  information  -  destroying,  nonlinear  recti¬ 
fier-integrator  (envelope  detector)  as  a  necessary  component  of  initial  processing.  The 
receiver  can  thus  operate  upon  the  spectral  phase  of  a  signal,  as  well  as  its  spectral 
amplitude.  In  terms  of  the  theory  in  Sections  2  through  4,  the  phase  as  well  as  the  mag¬ 
nitude  of  spectral  Taylor  series  coefficients  can  be  determined.  In  terms  of  matched 

56 

filter  theory,  true  pulse  compression  becomes  possible. 

A  disconcerting  property  of  Fletcher's  model  is  its  behavior  when  it  is 
utilized  as  a  sonar  signal  processor.  The  lov.  pass  filter  or  integrator  at  the  output  of 
each  critical  bandwidth  filter  results  in  an  output  pulse  that  is  at  least  0.  2  seconds  in 
duration,  regardless  of  the  bandwidth  of  the  input  signal.  The  range  resolutuu  of  a 
sonar  processor  equipped  with  such  an  envelope  detector  is  on  the  order  of  220  feet 
in  air,  or  900  feet  in  water. 

It  is  obvious  that  animal  echolocatior.  systems  have  much  better  resolu- 
lion,  in  fact,  J.  Simmons'  experiments  °  show  that  the  range  resolution  of  the  FM  bat 
fiptesicus  fuscus  is  predicted  by  the  autocorrelation  function  of  the  bat's  sonar  pulse. 
This  resolution  is  approximately  determined  bv  the  reciprocal  of  the  bat's  signal  band¬ 
width,  and  is  equivalent  to  about  half  an  inch. 
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8.3  —Continued. 


If  animal  sonar  data  is  to  be  explained  by  Fletcher's  model,  the  envelope 
detectors  must  i>e  altered  to  shorten  the  integration  time.  A  shorter  integration  time, 
however,  does  not  provide  an  adequate  explanation  for  the  measured  tone  duration  de¬ 
pendence  of  the  auditory  threshold. 

The  linear  model  does  not  suffer  from  the  above  problem.  Although  the 
impulse  response  of  each  5  Hz  filter  is  0.  2  seconds  long,  a  linear  combination  of  these 
responses  can  produce  a  pulse  with  a  much  shorter  duration.  The  minimum  width  of 
such  a  pulse  is  approximately  the  reciprocal  of  the  total  bandwidth  of  all  the  signals  that 
are  added  together. 

Envelope  detecting  the  output  of  each  5  Hz  filter  is  only  legitimate  if  one 
assumes  that  the  signal  is  a  tone  with  a  bandwidth  that  is  less  than  5  Hz  wide. 

The  new  model  allows  for  linear  signal  analysis  and  reasonable  range 
resolution,  but  docs  not  offer  a  completely  satisfying  explanation  of  critical 
bandwidth  data.  It  would  seem  that  a  minimum  number  (twenty)  of  the  N(f0)  combined 
5  Hz  filter  outputs  must  be  further  combined  when  narrowband  masking  noise  is 
introduced.  Perhaps  a  stimulus  with  a  bandwidth  exceeding  some  minimum  level 
elicits  power  series  spectral  analysis,  while  a  narrowband  tone  is  subjected  to 
ordinary  spectrum  analysis.  This  dichotomy  could  explain  the  divisions  of  many 
bat  signals  into  narrowband  tones  and  wideband  chirps. 


M  : 


ESL-FR115 


S. 


MINIMIZING  THE  VOLUME  OF  THE  WIDEBAND  AMBIGUITY  FUNCTION. 


The  volume  of  the  wideband  ambiguity  function  is  not  fixed  as  in  the  narrow- 
77 

band  case,  even  though  signals  are  defined  to  have  unit  energy.  It  is  still  possible,  how¬ 
ever,  to  define  a  set  01  signals  such  that  all  members  of  the  set  have  the  same  wideband 

57 

ambiguity  volume.  A  relevant  example  is  the  set  of  signals  with  Fourier  transforms 
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where  o  («;}  can  be  any  continuous  function  and  k  can  be  any  constant  greater  than  one. 
All  signals  whose  spectra  can  be  described  by  (S-l)  have  a  wideband  ambiguity  function 

O  f 

with  the  same  volume  as  the  narrowband  ambiguity  function  ,  i.e. , 
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where  E  is  the  energy  of  the  signal. 


For  good  range  and  doppler  resolution,  the  volume  under  the  ambiguity 
function  should  be  made  as  small  as  possible.  The  \  Glume  depends  upon  1  U  {u.*)  !  = 
A  (jjI,  since  ’’ 


Volume 


d  «.* 


(9-3) 


The  dependence  of  volume  upon  '  U  («.')  explains  why  o  (-)  can  remain  unspecified 
in  {9- 1)  . 
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An  undesirable  wav  to  minimize  the  volume  is  indicated  bv  the  ir.eaualitv 


58 


•*s2 


f  1  \ulT’  s)  |2  dtds  1  ^7?  TfJ 


1/2 


-.1/2  1/2 
X  (  s2  -  Sj  ) 


(6-4) 


where 


'-■L 


•  w  * 

[  A  <w)  ]2dw  i  ^  r  I  u  (o>)  ! 2  d  W  = 

O  0  -s 


03 

r< 


U(t)  !  “  d  t 


=  D 


(9-5) 


the  mean  square  duration.  If  Dt“  is  made  small,  the  volume  will  be  made  small.  For 

{  o 


a  gnen  energy,  however,  making  D.  small  implies  a  large  maximum  power,  a  quantity 
which  is  always  constrained  in  any  practical  system. 


A  meaningful  volume  minimization  problem  is  then  to  find  the  function 
/*<=  2 

A  (>u)  that  minimizes/  [  A  (u?)  /  w  ]  d  uj  as  in  (9-3),  with  constraints  that  (i)  E  . 

Jo  /•*  ■  _o  2 

the  signal  energy  be  kept  large,  and  (is) I  r  A  (m)  1  d  ,  a  lower  bound  for  D.”  . 
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The  functional  to  be  minimized  is  therefore 
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(9-6) 


where  A  and  \  are  Lagrange  multipliers  for  the  energy  and  time  duration  constraints, 
t  I 

The  Euler- Lagrange  equation  provides  a  necessary  condition  for  the  func¬ 
tion  A  (ce)  to  be  an  extremaloid  of  the  functional  (9-6)  : 
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9. _  _ --Continued. 


The  Legendre  necessary  condition  requires  that 
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(9-  9} 


in  order  that  the  extremaloid  he  such  that  the  functional  (9-6)  is  minimized.  Equation 
(9-9)  requires  that 


(9-1G) 


Equation  (9-8)  ?an  be  written  in  the  same  form  as  the  Coulomb  wave  equation 
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I  -( 2  77  /u?)  -  L(L*l)/ui  ]  A  -  0 


(9-10) 


To  make  (9-8)  identical  with  (9-10),  let 


X 


T 


1  /  2t?  : 


(9-11) 


The  two  equations  are  Mien  identical  if  L  -  0  or  L  -  1  .  The  solutions  are  the  Coulomb 
59 

wave  functions 


F  (??,«;)  and  F  .  t  77  ,  )  . 

o  -  i 
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Since  both  (r) ,  w}  and  F  ?  (77 ,  »*.')  are  solutions  of  the  same  differential  equation, 
they  must  be  equal,  indeed,  if  L  ~  0  is  substituted  into  the  recurrence  relation"^ 

L  [a  Fl  (  n  ,  )  f  d  x«]  -  (L  +  Ti  )1#  “  {??,  -?)  -  ( 7 j*  I.2/a*l  F  ( j)  ,  «.•)  , 

Li 

(9-12) 

one  comes  to  the  conclusion  that 


F0  (17  .  «■’ )  =  F  (  7? ,  w  ) 


(9-13) 


The  functions  A  (u.)  that  minimize  ambiguity  volume  while  constraining 
energy  and  time  duration  to  remain  large  are  therefore  Coulomb  wa-.e  functions  with 
L  =  0  . 

Since  A  (u.*)  is  defined  as  the  magnitude  of  the  spectrum  U  (u?)  ,  we  require 
that  A  (*)  c  0  .  Furthermore,  we  require  that  A  {->}  be  continuous,  so  that  the  integral 
in  (9-3)  is  defined.  These  two  implicit  constraints  imply  that 
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{9-14} 


where  z}  is  the  first  zero  of  Fq  ( 7] ,  )  for  w  >  0  .  For  example,  Figure  9-1  shows 
A  f— 1  for  1?  1  ,  ?7  =  5  ,  and  tj  -  10  .  For  energy  normalized  signals,  it  is  clear 
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9.  — Continued. 


from  (9-3)  that  the  volume  decreases  as  77  is  increased.  The  minimum  attainable 
volume  is  therefore  determined  by  the  maximum  frequency  that  can  be  obtained  with 
a  given  system. 


Equation  (9-14)  defines  the  amplitude  of  the  signal’s  spectrum, 
about  the  associated  phase?  From  the  viewpoint  of  volume  minimization,  th 


but  what 
e  phase  is 


arbitrary.  From  ;he  viewpoint  of  pulse  compression  techniques,  however,  a  nonlinear 
(chirped)  phase  function  is  desirable.  It  can  be  shown”  that  range- doppler  coupling 


will  be  avoided  if  a  log  phase  function  is  used  with  the  positive  semidefinite  functic 

* 

described  by  (9-14)  , 


*  to  general,  it  Is  shown  In  reference  5?  that  any  positive  semidefinite.  Analytic,  real 


frequency  function  ear.  be  multiplied  by  a  log  pat 


ion,  and 


resulting  transform 


will  correspond  to  a  signal  that 
in  Section  6,  where  some  Mvoti. 


=  no  range  doppler  coupling.  Another  example  Is  found 
gnals  were  matched  to  real,  positive  semidefinite 


spectra  with  log  phase  functicsis.  The  My  oils  signals  therefore  have  no  range- doppler 
coupling,  a  result  that  was  determined  experimentally  by  p_  A.  Johnson. 
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FURTHER  ANALYSIS  INVOLVING  THE  EQUIVALENCE  OF  AMPUTUDI 
MODULATION*  AND  LINEAR  FILTERING. 


Characterization  of  Time  Varying  linear  Svsti 


Is  Section  0 .  it  was  shows  that  amplitude  modulation  of  the  mused  signals 


u  (ll.  with  Fourier  transforms  gives  by  could 


~*e  vicweo  as  a  uncar  tillering  sroce: 


Heuristic  justification  for  this  equivalence  was  presented  is  Section  7. 


Figure  10- 1  si 


linear  system  which  can  be  made  equi*  aiest  to 


otaer  ttnear  system  over  a  given  irequency  interval.  The  differenitav^-rs  ':r?-  preceded  bv 
gams  that  determine  tbs  power  series  spectral  coefficients  of  F  (w)  ,  the-  svstem  transfer 
function.  To  characterise  a  time  varying  linear  system,  the  power  series  spectral  ec*?~ 
flcier.ts  are  allowed  to  change  with  time,  as  indicated  In  Figure  10- 1. 

When  3  mils c-d  signal  u  (t)  with  Fourier  transform  i2~8*j  is  introduced  at  the 
input  ot  me  s> stem  in  Figure  1C*-1 ,  me  modulation  functions  a  (tj  have  the  effect  shown 

S' 

m  Section  u.  In  other  wotds,  the  product  a  (t)  u  (ij  can  bt  described  as  a  linear  time 

n  ~ ■ 

Invariant  transformation  of  a{t). 

A  linear  time  invariant  irsasiontiation  can  again  be  characterized  by  a 
sequence  of  differentiators,  provided  an  appropriate  time  invariant  gain  is  associated 


iith  each  differentiation.  For  as  input  pulse  ileser?* 


1 2- St,  the  time  varying  linear 


l  ms  system  can  again  r>e  enamel- ri a*; d  oy  use  methods  in  Section  4.  Specif! cal!-. ,  oc 
can  design  a  time  invariant  receiver  to  detect  echoes  from  a  target  v  fees?  sla.ee  and* 
aspect  is  time  varying. 


tore  preo 


cisciv.  anv 


iy  otta  r  linear  system  with  a  transfer  function  ttet  is  analvtic  at  _  ~0. 


ne  argument  can  te  geri  railzed  to  otter  systems  by  including  integrators 


differentiators,  in  rigu 
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10.2  —Continued. 


One  way  to  transmh  binary  data  with  the  signals  (2-8)  is  to  use  linear 
increasing  period  modulation  for  one  signal  and  linear  decreasing  period  modulation  for 
the  other.  These  modulations  are  induced  by  using  a  frequency  phase  function  equal  to 
exp  [  ->  27m  log  w  /  log:  k  ]  and  exp  [  +  j  2m  log  u  /  log  k  ]  ,  respectively.  These  sig¬ 
nals  are  excellent  probes  to  determine  the  effective  transfer  function  of  the  channel  qua 
linear  system.  Therefore,  if  the  receivers  in  Section  4  are  used  (rather  than  only  two 
matched  filterst ,  one  can  characterize  the  channel  at  the  same  time  as  data  is  being 
trar»s  mitted. 


If  the  channel  is  stochastically  time  varying,  one  can  estimate  the  ex¬ 
pected  values  of  the  coefficients  for  a  linear  filter  model.  These  expected  values  are 

related  to  the  coefficients  A  that  describe  the  time  varying  gains  in  Figure  10-1,  i.e. , 
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(10-1) 


In  this  case,  a  (t)  is  a  stochastic  process  which  we  represent  by  the  expected  values  of 
•  n 

A  .  where 
nrn 
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(10-2) 


The  expected  values  of  A  are  related  to  moments  of  the  power  spectrum  of  i’ne  process 

nm  60 

a  (t).  For  a  stationary  stc  chaste  process, 

n  '  '  v 
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where  R  (r)  « — »  S  (u)  are  the  autocorrelation  function  and  power  spectrum  of  the 
an  an 

stochastic  process  a  (t), 
n 


10.  d _ Simplification  of  Complicated  Sonar  and  Communication  Problems  ;  A 

Return  to  the  Planar  Target/  Perfect  Channel  Situation. 

Radar/ sonar  operation  for  perfect  all- pass  channels  and  planar  ("point") 
targets  is  well  understood.  This  report  has  considered  nonplanar  targets  whose  shapes 
change  rapidly  with  time,  and  time  varying,  imperfect  channel  transfer  functions.  To  the 
extent  that  target  and  channel  transfer  functions  are  predictable,  appropriate  signals  and 
filtc-s  for  detection  of  a  target  and  analysis  of  its  transfer  function  have  been  derived. 

It  is  evident,  however,  that  there  is  considerable  unpredictability  in  detec¬ 
tion  of  a  target  from  a  specific  aspect  angle  at  any  given  time.  Stochastically  time 
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10.:?  — Continued. 


varying  targets  can  be  treated  via  the  analysis  in  Section  10.  2,  and  random  aspect  angles 
can  be  dealt  with  as  in  Section  4-8.  These  approaches  model  the  receiver  in  terms  of 
the  expected  values  of  the  random  parameters.  Although  performance  of  the  resulting 
system  may  be  optimized  in  terms  of  many  different  observations,  performance  may¬ 
be  poor  for  a  particular  observation,  e.  g.  ,  an  aspect  angle  such  that  he  observed  target 
parameters  are  very  different  from  their  mean  values.  In  terms  of  detection  theory,  the 
variance  of  the  sufficient  statistic  may  be  so  large  that  the  probability  of  false  alarm  is 
considerable  when  the  threshold  is  set  to  give  a  reasonable  probability  of  detection. 

Given  the  above  considerations,  the  sonar  theorist  (and  echoloeating  animals) 
may  yearn  for  the  simplicity  and  demonstrated  performance  obtainable  in  the  planar  target/ 
perfect  channel  situation.  One  way  to  achieve  this  simplicity  is  to  use  narrowband  signals. 
These  signals,  however,  are  of  little  use  for  range  resolution  and  suppression  of  clutter 
via  range  gating.  One  must  therefore  ask  the  question: 

"Does  there  exist  a  wideband  signal  such  that  the  response  of  its 
matched  filter  is  impervious  to  time  vary  ing  linear  transformations 
of  the  signal?" 


A  solution  to  the  above  problem  is  implicit  in  the  results  that  have  already- 
been  discussed  in  this  report.  In  Section  10. 1  it  was  demonstrated  that  time  varying 
transfer  functions  have  equivalent  time  invariant  counterparts,  if  the  functions  V  (w) 
in  Equation  (2-8)  are  used.  For  a  single  pulse,  any  time  varying  nonpianar  target  or 
channel  transfer  function  can  therefore  be  modelled  as  a  iime  invariant  linear  filter. 
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Furthermore,  it  was  shown  in  Section  2,  2  and  in  the  Appendix  tnat  if  the  parameter  k  is 

* 

sufficie,  tly  large,  i.e.  ,  if  lT  (w')  is  sufficiently  wideband  .  then  the  additional  echo 

o 

components  U  (u.)  -  al}  U  Its))  _  caused  by  time  invariant  filtering  are  uncor¬ 

related  with  C  (w).  lor  the  signals  defined  V  Equation  (2-8),  the  added  echo  components 
introduced  by  nonpianar  targets  and  imperfect  channels  do  not  affect  the  response  of  a 
filter  that  is  matched  to  U  {*•)  ,  if  k  is  large.  This  filter  is  matched  to  the  transmitted 
signal,  i.  e. ,  to  the  echo  that  would  be  convened  through  a  perfect  all-pas  >  channel  from 
a  planar  target. 


m  i- 
it 


1  I 


For  the  Rivals  (2- Si  with  sufficiently  large  k,  the  response  of  a  single 
matched  filter  is  therefore  impervious  to  time  van  ing  linear  transformations  of  the 
signal.  Echoes  can  therefore  be  processed  and  interpreted  in  the  same  simp’e.  straight¬ 
forward  way  that  would  apply  to  planai  "arget/  perfect  channel  situations. 


For  communication  systems,  utilization  of  the  signals  (2-8)  should  result 
in  minimum  degradation  of  matched  filter  performance,  e\en  if  the  channel  is  rapidh 
time  \ arcing  and  unequalized.  The  wide  bandwidth  allows  for  resolution  of  multipath 
receptions. 


|I 
1  I 

i  i 


*  Note  that  the  isolation  between  echo  components  increases  with  bandwidth.  For  ihc- 
signals  in  Equation  (2-8),  the  planar  target  model  becomes  more  realistic  as  bandwidth 
is  increased,  rather  than  decreased. 

40 

**  This  result  explains  an  observation  made  by  E.  I,  Titlebaum  and  R,  A.  Johnson*  , 
\iz. .  that  a  Mvotis  pulse  is  nearly  impervious  to  low  pass  filtering,  so  far  as  matched 
fiber  oerformanee  is  concerned. 
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10.  !  Application  to  Diagnostic  Ultrasound. 


If  the  wrong  signals  and  filters  are  used,  it  is  easy  to  see  how  a  wider  pulse 
bandwidth  could  introduce  more  confusion  than  it  eliminates,  especially  if  the  planar 
target/ perfect  channel  condition  is  explicitly  or  implicitly  assumed. 


An  example  cf  this  confusion  is  to  be  found  in  diagnostic  ultrasound  research* 
When  a  constant  f requeue},  pulse  (the  fa  'orite  signal  for  diagnostic  ultrasound  ’  ’  ') 

is  made  more  narrow  so  that  its  bandwidth  increases,  "any  difference  in  amplitude  or 
velocitv  for  the  different  frequency  components  within  this  bandwidth  will  distort  the 
*  filtered  :  pulse  and  generally  widen  it.  "  5  ’ 


Diagnostic  ultrasound  techniques  are  generallv  used  to  construct  an  image 
of  an  object  by  means  of  multiple  reflections  from  the  reflector's  surface.  For  most 
wideband  signals,  this  image  is  degraded  b\  imperfect,  sound  channels  (intervening  tissue) 
or  b\  e:  nature  of  the  reflecting  surface.  Tiiis  degradation  occurs  because  of  distorted 
pulses  at  the  .Hatched  filter  output.  The  distortion  can  be  manifested  as  widening  of  the 
output  pulse,  unwanted  ambiguity  peaks,  and  a  displaced  maximum. 


It  has  been  shown  (Section  10.  3)  that  the  increased  resolution  inherent  in 
wideband  signals  (for  the  planar  target.' perfect  channel  situation)  can  still  lie  obtained 
under  imperfect  conditions.  A  wideband  signal  can  be  used  for  ultrasonic  diagnosis  if 
the  corresponding  matched  filter  response  is  impervious  to  linear  transformations  of 
the  signal.  Equation  (2--?)  describes  the  Fourier  transform  of  such  a  signal. 


*  The  dispersive  effect  of  different  propagation  velocities  at  different  frequencies  can  be 
modelled  as  a  pulse  compression  or  expansion.  The  Doppler  tolerance  of  the  signals 
(2-m  and  their  lack  of  range- Doppler  coupling  (Section  9.  last  paragraph*  ensure  that 
that  such  effects  'will  not  dvgrndo  the  output  of  a  filter  that  is  matched  to  the  transmitted 
wjixpfnrm. 
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10.  5 _ The  FM  Equivalent  to  any  Measured  Transfer  Function  *,  FM 

Demodulation. 


In  Section  G,  amplitude  modulation  functions  were  written  as  power  series 

with  complex  coefficients.  A  frequency  modulation  function  can  be  written  in  the  same 

way.  It  would  therefore  seein  that  the  receiver  in  Figure  6-7  can  be  adapted  to  estimate 

/ 

the  coeffic.^nts  of  a  function  0  (t)  that  is  used  to  frequency  modulate  the  signal  u  (t). 


The  frequency  modulation  process  multiplies  u  (t).  the  inverse  Fourier 
transform  of  one  of  the  functions  in  (2-S),  by  e.xp  [  i  o(t)  j.  Although  multiplication  of 
u  (t)  by  exp  [  j  0  (t)  jean  be  interpreted  as  an  amplitude  modulation  process,  a  complex 
representation  of  u  (t)  gives 


i  0  (t) 


u  (t)  =  e 


jp(t) 


[  I  u  Cl) 


jo  (t) 
e  u 


u  (t)  i  e 


j  l 


OuW 


o  (t)  , 


(10-4) 


i  he  phase  function  associated  with  u  (tl  has  an  added  information-carry- 

/ 

ing  term.  For  standard  FM,  (d/dl)  [o  (t)  +  d  (t)  ]  -  u  +  6  (t)  ,  where  cc  is  the 

t  U  O'  o 

carrier  frequency  and  d  it)  is  the  modulating  signal. 


The  desired  information  is  contained  in  the  power  series  coefficients  of 

c  (ti  .  The  demodulation  process  should  therefore  provide  estimates  of  D  ,  D  ,  D  , 

o  1  2 

...  ,  D  ,  where 
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10.)  — Continued. 


In  Section  G.8,  a  method  of  finding  the  coefficients  A  corresponding  to  an 
amplitude  modulating  function  was  given.  In  terms  of  a  modulating  function  exp  [j  e(t)  j, 
we  have 


c  i  e  t‘) 


X 

n-Q 


n 


A  t 


n 


(10-6) 


The  problem  is  to  find  a  transformation  that  relates  the  .A  coefficients  in 
(10-6)  io  the  I)  coefficients  in  (10-5).  Differentiating  (10-6)  gives 


±  r  cj  e  <*>  1 
dt  L  - 


X 

n=r 


A  t 
n 


n-1 


(10-7) 


where 


B  -  (n  +  1)  A 
n  n  +  1 


(10- S) 


In  matrix  notation,  (10-8)  can  be  written 


B  =  H  A 


(10-9)  (10-9) 


where 
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Evaluating  the  left  s»ds  of  (10-7), 
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10.  ;•  — Continued. 


(10-15)  is  implemented.  Although  the  demodulation  is  accomplished  by  a  matrix 
transformation,  the  transformation  is  nonlinear  because  is  processed  with  a  matrix 
whose  elements  depend  upon  A. 

Notice  that  the  inverse  relation,  A  =  j  [  T  ^  1  H  „*  5  D  ,  docs  not  make 

'  _1 

sense,  since  A  must  already  be  known  in  order  to  calculate  T  .  This  observation 
suggests  that  the  demodulation  process  to  obtain  the  coefficients  of  6  '  Jit  is  irrevers¬ 
ible  (information  is  destroyed).  Given  D,  one  cannot  uniquely  reconstruct  the  A  coef¬ 
ficients  which  correspond  to  D. 

Equation  (10-15)  demonstrates  that  a  linear  filtering  operation  can  be 
interpreted  as  a  frequency  modulation  of  the  signal  u  (t)  ,  provioea  the  right  hand  side 
of  (10-15)  exists. 
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1 0.  S  Taylor  Expansions  about  Frequencies  Other  than  2ero. 

It  has  beer  shown  that  amplitude  modulation  and  linear  filtering  can  have 

identical  effects  upon  spectral  Taylor  coefficients.  These  coefficients  can  be  defined 

in  terms  of  a  Tavlor  expansion  about  an\  arbitrart  frequency  ^  ,  as  well  as  about 

o 

_  o.  Returning  to  Section  4.  4,  we  can  redefine  the  functions  U  {--)  and  F  (*•)  in  terms 
of  expansions  about  j  : 


V  {s) 


N 

I 

n-D 


u  P 

n 


F  (cc)  = 


1 

n--0 


P  - 


a’  ) 

o 


I 

n  -0 


N 

2 

n~0 


U  wt* 

n 


n 


Now  if 


N 


B  p)  -  F  p)  U  p) 


2  *. 


n  0 


~  n 

>  b  p  -  } 

Z,  n  ° 


n=o 


it  follows  that 


b  i  f 


(10-16) 


(10-17) 


HO- IS) 


(10-19) 


10-15 


-Continued. 


where  T  is  the  suuare  matrix  in  (4-33). 
-  u  • 


Furthermore, 


b  *■  T- 


(10-2 


so  that  the  same  resul:  applies  to  coefficients  of  a  Tavlor  expansion  about  a>  £  0 

o 


Similarly,  if 


a  (t)  = 


2 


(10-2 


as  in  (6- SI,  then 


a  (t)  u  (t) 


(10-2 


has  coefficient 


I  ? 


where  c  is  proportional  to  the  n—  derivative  of  C  (x)  evaluated  at  u.-  -  •»:  .  In  othe 
n  o 

words,  (6- 11)  becomes 


Channel  transfer  functions  F  (oj)  that  satisfy  the  above  inequality  will  not 
seriously  affect  the  matched  filter  response.  Notice,  however,  that  Titlebaum’s  counter 
example  does  not  satisfy  the  inequalitt’.  Condition  (10-30)  can  therefore  he  used  as  a 
test  to  determine  the  applicability  of  the  theory  to  any  given  situation. 
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11. _ SYNTHESIS,  TRANSMISSION,  AND  UNDERSTANDING  OF  SPEECH. 

Applications  of  the  foregoing  ideas  to  speech  processing  became  obvious 
in  the  light  of  two  observations: 

1.  The  human  and  auditory  spectrum  analyzer  is  apparently 
set  up  to  perform  the  operations  discussed  in  Sections  3 
and  4.  Specifically,  signals  are  analyzed  to  determine 
their  spectral  Taylor  series  coefficients.  Receiver  oper¬ 
ations  are  greatly  simplified  when  (i)  the  coefficients  of 
interest  are  associated  with  a  filter  transfer  function  F(cu), 
and  (ii)  the  received  signal  has  Fourier  transform  F  C&u)  x 
U  (w),  where  U  («)  is  defined  as  in  (2-8). 

2.  Human  speech  can  be  modeled  as  the  passage  of  a  pulse 
train  or  white  noise  through  a  slowly  time  varying  linear 
filter,  so  that  the  information  content  of  speech  is  obtained 
through  analysis  of  the  filter  transfer  function. 

6T  6^- 

The-  usual  speech  synthesizer  model  **  '  is  shown  in  Figure  11-1.  An 
equivalent  model  is  shown  in  Figure  11-2.  In  Figure  11-2,  the  white  noise  generator 
has  been  replaced  by  a  generator  of  pulses  that  occur  at  random  times  (Poisson  process 
or  shot  noise)  and  each  fillet  (for  voiced  and  unvoiced  components)  has  been  divided 
into  a  filter  with  transfer  function  U  (w)  and  a  second  filter  in  cascade.  Knowledgi  of 
the  transfer  functions  of  the  second  filters  (F  (u:)  and  F  (u)  in  Figure  11-2)  implies 

V  u 

knowledge  of  Lhe  information  content  of  the  sounds  produced  by  the  speech  synthesizer. 


Generator  of 
a  Train  of 
Impulses, 

T  Seconds 
Apart 


Filter  for 
Voiced 
Component 


SL-PRUd 


Fixed 

Filter 


White  Noise 
Generator 


Filter  for 
Unvoiced 
Component 


Figure  11-1.  A  Model  for  Speech  Synthesis,  Devised  and  Tested  by  Flanagan, 

et.  al.61 
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Figure  11-2, 

A  Model  for  Speech  Synthesis  that  is  Equivalent  to  the  .Model 
Figure  11-1. 

U  (a?)  is  given  by  Equation  (2-8). 
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It  was  deduced  in  Section  3  that  the  ear  is  designed  in  such  a  wav  as  to 

characterize  tie  transfer  functions  of  the  filters  F  (u?)  and  F  (w)  (Figure  11-2)  in  terms 

v  u  * 

of  their  power  series  spectral  coefficients..  Furthermore,  a  specific  method  has  been 
dented  for  processing  the  data  at  the  output  of  the  ear’s  critical  bandwidth  filter  bank. 
Receiver  configurations  to  detect  a  given  set  of  coefficients  (by  implementing  a  likeli¬ 
hood  ratio  test)  or  to  estimate  the  coefficients  of  an  unspecified  filter  fbj  maximum 
likelihood  estimation)  wore  deduced  in  Section  4.  A  complete  receiver  specification 
has  thus  been  obtained  For  both  detection  and  estimation. 

The  above  deductions  could  have  been  made  without  the  added  insight  that 
has  beer,  provided  by  animal  c-eholocaiion  studies.  Nevertheless,  It  has  been  observed 
that  signals  with  transform  U  (w)  (Equation  2-8  and  Figure  11-2)  are  actually  utilized 
bv  both  bats  and  cetaceans  to  characterize  a  linear  filter  (i.e. .  a  nonplanar  sonar 
target).  This  observation  was,  of  course,  the  impetus  for  an  interpretation  of  the 
auditory  spectrum  analyzer  as  a  device  for  power  series  spectral  characterization. 

More  Important,  animal  echolocation  signals  provide  evidence  that  the  speech  synihe- 
sizor  in  Figure  11-2  is  more  than  a  convenient  mathematical  model.  Animal  sonar 
strongly  suggests  that  the  human  brain  interprets  speech  exactly  as  if  phonemes  were 
synthesized  by  the  model  shown  in  Figure  11-2. 

Ammals  that  ucholoeate  are  forced  to  utilize  sonar  signals  that  are  com¬ 
mensurate  with  their  built-in  receivers.  The  same  reasoning  implies  that  humans 
are  also  forced  to  utilize  -.ommunication  signals  that  are  commensurate  with  their 
built- in  sound  analyzers.  The  optimum  processing  methods  of  Section  4,  together 
with  critical  bandwidth  illtei  transfer  functions,  should  therefore  yield  an  extremely 
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11.  _  —  Continued. 

simple  characterisation  o?  speech.  A  reasonable  hypothesis,  then,  is  that  if  Taylor 
series  spectral  analysis  is  applied  to  speech,  the  basic  information-earning  parameters 
of  the  speech  waveform  should  be  obtained.  Access  to  these  parameters  would  greatly 
simplify  the  computer  understanding  of  speech,  as  well  as  data  compression  for  trans¬ 
mission  of  speech. 

Finally,  one  can  apply  the  foregoing  ideas  to  a  very  interesting  topic  -  the 

analysis  of  animal  communication  signals.  By  determining  the  appropriate  coefficients 

of  cetacean  communication  sounds,  one  can  compare  the  data  rate  with  that  of  humans. 

The  analysis  may  also  provide  a  method  for  computer  understanding  of  cetacean  speech 

* 

(if  there  is  such  a  thing  ) ,  as  well  as  that  of  human  beings. 


*  Some  observations  concerning  exchange  of  information  between  cetaceans  are  discussed 
on  pages  55-58  of  R.  Ster.uit's  book. 
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12. _ SUMMARY. 

A  particular  set  of  signals  provides  an  excellent  match  to  the  measured 
echolocation  clicks  of  both  Tursiops  truncatus.  the  Atlantic  bottlenose  dolphin,  and 
Myotis  lucifugus,  the  little  brown  bat.  This  signal  set  can  be  derived  in  at  least  two 
different  ways.  The  derhations  are  based  upon  observed  properties  of  sona  targets, 
the  animals’  auditory  systems,  or  the  measured  signals  themselves. 

The  derivations  incorporate  new  ideas  concerning  sonar  target  characteri¬ 
zation  .  mammalian  auditory  processing,  and  modulation  theory.  Each  new  insight  is, 
in  turn,  pregnant  with  mans  applications  that  seem  far  removed  from  the  rather  esoteric 
desire  to  formulate  a  unified  theory  for  animal  echolocation. 

The  first  derivation  is  based  upon  the  concept  that  sonar  targets  can  be 
characterized  by  their  transfer  functions,  just  as  in  linear  system  theory.  An  impli¬ 
cation  of  this  concept  is  that  proper  signal  and  filter  design  can  accomplish  automatic 
clutter  rejection  and  target  recognition,  on  the  basis  of  reflector  shape.  A  separate 
theory  of  signal- filter  design  for  signal  to  interference  ratio  maximization  can  thus 
be  applied,  independent  of  animal  echolocation  systems. 

The  second  derivation  is  based  upon  the  concept  of  a  linear  filtering  system 
to  determine  the  spectral  Taylor  series  coefficients  of  any  input  signal.  Once  the  system 
has  been  mathematically  defined,  it  is  hard  to  ignore  ns  similarity  to  psyehoaeousiic 
models  of  the  mammalian  auditory  processor.  Furthermore,  the  system’s  operation 
(based  on  likelihood  ratio  testing  and  maximum  likelihood  estimation)  is  very  simple. 


provided  that  the  processed  signal  spectra  can  be  written  as  a  product  F  (_l  U 
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12.  — Continued. 

F  {w}  is  unknown  and  U  (w)  belongs  to  the  set  that  has  already  been  matched  to  bat  and 
porpoise  echolocation  eiieks.  Obviously,  if  U  (a?)  is  used  as  a  transmitted  sonar  pulse, 
all  echo  spectra  have  the  required  form. 

The  above  argument  not  only  constitutes  an  alternate  derivation  of  tat  and 
dolphin  sonar  signals;  it  also  provides  insight  into  the  workings  of  the  human  auditory 
system.  Particularly  important  is  the  assertion  that  human  speech  is  both  synthesized 
and  characterized  by  the  Taylor  series  coefficients  of  F  \w),  where  the  speech  wave¬ 
form  itself  can  be  written  as  the  product  F  (w)  U  (w). 


In  a  recent  article  in  National 


A.  Novick  stated  that  "research 


on  these  highly  specialized  mammals  f  bats  J  may  help  explain  how  man’s  brain  proces¬ 
ses  information.  The  basic  design  of  a  brain  -  be  it  man’s  or  bat’s  -  is  much  the  same. 
By  tracking  the  processing  of  sound  through  a  bat's  brain,  which  is  much  simpler  than 
man’s,  we  hope  to  understand  how  sound  is  coded,  analyzed,  integrated,  and  acted  upon." 
The  second  derivation  indicates  that  Dr.  Novick’s  hopes  tote  a  firm  foundation  In  mathe¬ 
matical  modeling  as  welt  as  in  physiology. 


Taylor  series  spectral  analysis  leads  to  another  interesting  result*  Multi¬ 
plication  of  two  time  functions  or  of  two  frequency  functions  can  to  described  by  ibe  same 
tvpe  of  transformation  in  the  frequency  domain,  if  Taylor  series  coefficients  are  used  to 
describe  the  sigrsals.  For  the  signals  used  by  tots  and  porpoises,  multiplication  by  a 
function  in  the  time  domain  (modulation)  and  multiplication  by  another  function  In  the 
frequency  domain  (linear  filtering)  car.  have  the  same  effect.  This  result  can  be  applied 
to  characterization  of  time  varying  linear  systems.  Such  an  application  Is  undoubtedly 
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useful  to  sonars  that  must  contend  with  targets  such  as  flying  Insects,  i.  e. ,  targets  with 
time  varying  transfer  functions.  It  is  not  surprising  that  a  system  which  treats  ampiitsd 
modulations  and  l.rr.e  Invariant  linear  transformations  In  the  same  way  should  be  useful 


lor  tne  am 


lysis  of  time  varying  linear  traasformatlcgis.  "Every  linear  Insiantsueoii 


system  may  he  characterised  by  a  relation  of  the  form  v  (1)  A  (t)  s  (t; 

71  " 

is  constant,  the  svst  em  Is  i  *xco. 


if  a  m 


The  above  characterization  lends  to  a  straightforward  way  of  equalizing  a 
randomly  time  varying  channel  with  a  linear  filter.  Still  another  application  Is  the  are 
of  matched  filtering  to  demodulate  an  amplitude  modulated  pulse. 

The  signal  set  which  seems  to  be  basic  to  animal  echclocaiica  is  apparently 
hash,  to  mar.y  aspects  of  enginee  ring.  The  new  signals  can  be  used  as  robes:  communl- 
cation  waveforms.  If  tiiey  have  Tallin  -lent  bandwidth,  the  signals  yield  correlator  outputs 
that  are  not  degraded  bv  time  varying,  jmoerfeci  chanr-els. 


An  interesting  \u 


of  the  new?  signals  is  that  thev  are  broadened  {their 


meats  square  bandwidth  is  increased)  when  they  are  multiplied  by  air  function  In  the 
freouenev  domain  that  is  snalvUc  at  ~  -  0  .  In  this  fr-nse.  the;  seem  to  vossess  a  ~ 


•-  urt  ti  o  ropes 


Equation  (2-Aj  constitute  the  most  emeient  prates  for  eha racte sizing  a  lir 


easure  ci  e  -i  is 


:icr,cv  Is  toe  miarmattoR 


is  aooesrea  so*  a  given  it 


sauare  signal  bandwidth. 


The  derivations  hinge  upon  Ire  existence  m  set  of  deeorrelated  signals  such  that  esc! 


,  HkVdl'sfW ’fl  **  “Ml IMWil^^4KillUI9ili4|ti^^iWiRl#:|HiH^]l;)l|JR^|fi^^l^i#Ij^HVW|UK|<MHi!4!^^  M4b  it**  »!!| 


u  Ml,  ll  V, 


difficult  to  establish  the  uniqueness  ot  ssea  a 


however,  that  decorrelation  is  a  more  stringent  condition  to  the  us 
concept.  Orthogonality  Implies  zero  inner  product  fur  Kro  iranslatf 


u  ortco-ssan 


cut  itmctiQ 


resrect  to  another.  Decor 


Marios  requires  verv  smau 


nr^*iic!  tor  sji 


ssble  translations.  At  the-  present  time,  the  author  Is  pleased  to  hare  found  only  one 
solution  to  the  problem.  Although  other  solutions  may  exist,  it  would  appear  that  the 
author’s  result  compares  nicely  with  the  sole- ions  that  are  implemented  is  the  mam* 
mall  an  auditory  processor  and  in  the  ecitolocatuffi  signals  u  sed  or  bats  arm  osipoBss, 

Related  topics  that  have  been  discussed  in  this  report  Include  a  new  monel 
for  the  auditory  processor,  the  volume  of  the  wideband  ambiguity  function,  and  clutter 
suppression  with  a  completely  linear  receiver. 

The  new  model  of  the  auditory  processor  was  motivated  by  the-  realization 
th^*  a  well  accepted  defer  mode:  cannot  explain  animal  evholocafion  performance.  The 
now  model  attempts  to  reconcile  tbs  concept  of  critical  bandsridth  spectral  analysis  with 
neural  nrocossing  and  behavioral  observations. 


An  interesting  result  concerning  ambiguity  volume  has 


been  cte  rived,  some 


of  the  signals  matched  to  animal  cch-olccation  waveforms  rave  wtoeomd  am^igai'.y  voiunu 


to  the  volume  under  the 


ambiguity  function.  Signals  that  minimize 


wideband  ambiguity  volume  under  realistic  constraints  (duration  and  era 


rgy  tvOui  iStg«j 


Iso  been  derived.  The  Importance  of  these  signals  to  animal  ecfeolocstion  theory 


has  vet  to  be  determined. 


.PHI' 


men. 


Clatter  sappresslqo  with  a  square-  !s«  detector  at  the  oslpai  -A  a  linear  filter 
leads  to  the  theory  reviewed  fa  Section  IS,  A  dtfier^t  approach  caa  be  applied  to  a 
receiver  that  does  not  use  square  lap  detection.  The  new  approach  mcHcsrs  iaai  aai 
bane  sinusoids  constitute  an  unacceptable  solution  to  the  clatter  sspresslon  problem,, 
unless  the  target  is  moving  with  respect  So  the  clutter.  A  mare  robust  solution  Is  deter 
mined  by  uslug  a  basis  set  of  wideband  signals.  Us  resulting  solution  is  less  affected 
bv  unforeseen  DsriyrbaiioBS  of  the  esviresment  or  the  sonar  emtipmeai. 


The  above  summary  Indicates  that  most  of  the  questions  ashed  is  Section 
1. 1  have  indeed  been  answered.  Special  importance  can  be  attached  to  the  following 
results; 


Sonar  Systems. 

linear  niter  theory  provides  a  straightforward  method  of  exploiting 
the  phenomenon  that  wideband  echoes  often  are  different  from  the 
traasm.rted  signal. 


Through  proner  design,  die  pnenomenon  can  oc  igttorea  or 
ci ream vanted  with  little  degradation  In  matched  filter  response 
(Section  10.  3i . 

Developme  nt  of  the  theory  leads  directly  to  system  implemeatailffiis 
to  achieve  automatic  target  recognition  and  clutter  rejection. 
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acoustic  osts,  ctsscesslBf  critical  basd^KSh, 

Tiss  theory  seems  Is  imply  the  esSsteae©  of  a  slightly  different  bb^I 
of  the  SHBtow  processor  tsm *  t m  cr,e  feat  is  currently  In  rogae,  Use 
results  also  stress  the  importance  sf  efferent  signals  f by  neural  proces¬ 
sing  of  SKiItory  caia.  These  observations  "stmM  seem  to  suggest  some 
important  sot-  experiments  In  searopinrslolsgf  sw  ssjofcosceestics. 


to  It  darsftriisw  :i  ?i— s  v:  rrieg  I tor  st^bs  Is  ih  r.r^ssr.; 
ef  ssise.  Otfcer  »sp3Ie«S«is,  to  sses  areas  as  wjoerlasu  rs^r. 
sstesseoic  gskiaace  for  fite  bliss!,  aas  sslssdc  geopferslcsl  prospeei- 
mgu  ssosM  be  strafgifiIbrsar*L  Fsriber  BgSc^_as  are  ielt  to 
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APPENDIX 


t im 


A.  1 _ Derivation  of  Signals. 

The  theory  in  this  report  is  based  upon  the  idea  that  an  echo  from  a  sonar 
target  css  be  expressed  as  a  linear  transformation  of  the  transmitted  signs!.  Targets 
can  therefore  be  characterized  by  their  transfer  functions,  just  as  linear  filters  are 
characterized  in  electrical  engineering.  The  target  transfer  function  can  be  parrmeiei- 
ized  by  expanding  the  function  in  a  Taylor  series,  i.  e. ,  a  polynomir.al  in  frequency 
whose  coefficients  describe  the  target.  One  is  then  faced  with  the  task  of  -determining 
the  Taylor  coefficients. 

In  order  to  estimate  the  coefficients  with  a  minimum  of  mathematical  data 
manipulation,  one  should  use  a  signal  that  becomes  uncorrelated  with  itself  when  its 
Fourier  transform  is  multiplied  by  w11 ,  n  =  0,  1,  ...  ,  N.  If  such  a  signal  exists,  a 
bank  of  filters  can  be  used  to  decompose  echoes  into  components  whose  amplitudes  are 
the  desired  target  coefficients. 

One  set  of  functions  that  become  uncorrelated  when  they  are  multiplied 
by  t?  ‘  has  the  property  that 


O  (w)  = 


u  (w/kni 


(A- 1) 


together  with  the  property  that 
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A.  1 


— Coniinuec 


*n. 


X™  <T-k  ) 


k 


B/ 


( 


U  (uj)  U*  (61/ k  )  e  ' 


o 


<<  i  X,,..  1}  I?,  for  ail  r 

Uu 


{A- 2) 


In  Equations  (A- 2)  and  (A-  2),  U  (w)  is  the  Fourier  transform  of  the 
.  n  1 2 

Analytic  signal  u  (i)  and  (  (F,  k  )  I  is  the  wideband  ambiguity  function  of  (he 

signal.  The  ambiguity  function  describes  the  response  of  the  filter  U*  (&>)  to  the  energy 
normalized  signal  k  U  (w/  k  *) . 


To  derive  the  s.g^als  that  satisfy  Equation  (A— 1) ,  we  consider  a  more 
general  condition; 

if  U  (fa?)  =  C  (P)  U  (t.-/kP)  (A- 3) 

where  p  is  any  real  mimber.  Differentiating  both  sides  of  (A- 3)  with  respect  to  P  , 
(log «)  eP  log "  U  («)  =  C'W  «;u/kP) 

-  C  CP)  W  (log  k)  e  P  Iog  K  Uf  (6J  e"  v  log  *S 

(A-4) 


A- 2 


ES] -~PP.il  o 


rsitlmied. 


Setting  P  =  0  in  (A-4)  gives 


U  ;  (fa» 

O  i>5) 


c  MO)  -  log  < 

C  (0)  — '  log  k 


(A- Si 


By  setting  P  =  0  in  (A-3)  ,  it  is  evident  that 


(0)  -  1  . 


(A-6) 


Integrating  both  sides  of  (A- 5)  and  using  (A-6) . 


,  C  ’  (0)  log  os  i  toe  w'  ~ 

log  U  (os)  -  -  •  V-"-  - 

log  k  2  log  k 


-  co.'stant 


U  (csj)  =  (constant)  w  e 


-  (log  -  ) 

2  log  k 


(A- 7) 


For  the  special  case  p  -  n  =  Integer,  U  («i)  can  be  multiplied  by  any  function  F  («.*)  w  i 
the  orone  rtv  that 


I V«  *  .|«» ill" MMMMM'J  IIWI— ii»W MMM HM  >, >M 


ii 


i  s 
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Ambiguity  Function. 


Tlie  wideband  ambiguity  function  of  an  Analytic  signal  with  Fourier  trans¬ 


form  U  (a;)  and  =  1  is 


X  (  T,  s  )  | 
mu  1 


1/2  f  U(a 

it  s  J 


_  i  T  9 

(CO)  U*  [Co/  s)  e  J  d  w  !  "  . 


{A- 13) 


For  the  functions  (A- 12)  , 


it  .  /,  n  ,  -TiV  -  n  /2  n 
U  (w/k  )  =  k  '  to 


U  (co) 


(A- 14) 


so  that 


w”  U  (a;)  =  knV  +  n/2  u^/k®) 


(A- 15) 


C  =  k 
n 


ny  +  n  /2 


(A- 16) 


in  (A-l). 


Using  {A- 14)  , 


-2m  v  -  2m  -m 


V  (T  k2m)  I2  -  1  k  1  f  2m  2  -jcoT  ,2 

Auu  ’  k  !  ~  !  - ”*  I  w  i  U  (w)  !  e  d  u.'  \ 


(A- 17) 
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A,  2  — Continued. 


Therefore. 


max 

_co  <£  j-  ro 


x  (T.k2m)  r 

^uu  v  ' 


=  k 


■2m 


(A- 21) 


In  Equation  (A- 21)  ,  max  |  X  (T»  k^m )  }  "  is  the  maximum  output 

T  2m 

power  of  a  filter  with  a  scale  mismatch  of  k  .  For  example,  the  signal  could  have 
Fourier  transform  k  ^  2  U  (c o/  kD)  ,  while  the  transfer  function  of  the  filter  could  he 
matched  to  a  different  signal  with  scale  factor  kn  +  2m,  i.  e. ,  k~  ^  y  (oj/ kn  + 
Equation  (A- 21)  gives  the  maximum  power  at  the  output  of  the  filter  when  the  mismatched 
signal  is  applied  to  the  input. 
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Mean  Square  Bandwidth. 


The  mean  square  bandwidth  (MSBW)  of  a  function  is  defined  as  the  difference 


(MSBW)  =  D/; 


2 

«  . 


where 


CD 

o>2  !  U  (u)  !2  d«  /  Eu  (A- 22} 


CO 

/  " 


|  u  (o>)  |  d  U!  /  E 


(A- 23) 


S3 

*/ 


{  U  (w)  |  dw 


{A-  24) 


From  (A- 15)  , 


o  9  9  9 1/  +  1  . 

«2  i  U  (u)  =  \ta  U(«)l  =  k  !u(«/k)  | 


(A- 25) 


U(«)  I2  =  |«1/2U(«)  i2  =  kV+1/4  |U(«/k^  i2 


(A- 26) 
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A*  3  — Continued, 


In  {A- 26),  it  is  assumed  that  the  periodic  function  G  (log  w  /  log  k)  in  (A- 12)  has  the 
property  that  G  (x  -  n/  2)  =  *  G  (x) .  Using  (A- 25)  and  (A- 26)  to  evaluate  (A- 22)  and 
(A- 23) , 


’M 


CD 

'/ 


j  U  (o?/k)  I  dtu  /  E 


u 


2^  -f  2 


(A- 27) 


1/11/4  J  J  U  (cVkl/2)  | 2  d  u;  /  E 


k^3/4 


(A- 28) 


Therefore, 


(MSBW)2  ,  k2""2-  k2^3/2  =  k2**372  (kl/2.  „  .  (A- 29) 


For  v  =  0  , 


MSBW  =» 


k3/4  <k1/2-  ),1/2 
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Approximate  Fourier  Transforms. 


It  is  difficult  to  obtain  exact,  closed  form  expressions  for  the  inverse 


Fourier  Transforms  of  the  functions  {A- 12),  where,  for  example, 


G  (log  w  /  log  k)  =  e 


-j2  «  n  log  a)  /  log  k 


(A-  30) 


It  is  relatively  easy,  however,  to  obtain  approximate  expressions  for  u  <t)  by  using 

T>  73 

the  stationary  phase  principle.  “* 

First,  we  obtain  an  integral  with  limits  (-®  ,  00  >  by  changing  variables  ; 


u{t)  =  T7~  I  " 


v  -  (log  w)  /  2  log  k  -  j2  »  n  log  u/ log  k  icet 
e  e  e  aw 


j  r®  (v  1)  -w^/  2  1ogk  -j2  “n  Wj/  log  k  ite^l 


*7. 


e  do 


(A- 31) 


According  to  Papoulis  (ref.  73,  pp.  140  -  141)  , 


U  to)  ejtfJ  {U:)  dcu 


U  toQ)  j[tptoQ)  *  Sgn  p  'toQ)j 


fe*  t  |  fJ.  {Ci?  ) 


(A-  3  2) 


Sgn  ft  (u  )  = 
\) 


*1  ,  u  '  too)  s  0 
-1  ,  <0 


{A- 33) 
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A.  4  — Continued. 


and  where 

w  =  the  stationary  point  of  p  (o>),  i.  e. ,  the  value  of  us  where  u  ‘  (a?  )  =  0. 
o  o 


For  tt»e  expression  in  (A- 31) 


P  {us) 


2?ngj 

tlogk 


us 


-  Z  US  /  t 


(A- 34) 


where 


2  ~  n 
2  =  logk 


Then 


u) 

II  '  {us  )  =  e  0  -  z  /  t  =  0 
o 


i  A- 35) 


so  that 


us  =  log  {z  /  t) 
o 


(A- 36) 


and 


us 


p*(oj)  =  e  *  z/t  =  |  at/  t  1  Sgnt.  (A-37) 


A-li 
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A*  4 


— Continued 


From  (A-  37),  we  see  that 


Sgn  \i '  (w  )  =  Sgn  t  if  z  >  0. 


Substituting  (A-34)  -  (A- 33)  into  (A- -3 2),  we  obtain 


u  (t) 


(v  4  1)  log  (2«  n  /  t  log  k)  -  [log  {2*n  /  t  log  k)j  /  2  log  k 

|  23-n/ tlogk  | 


.  e 


n 

jt  [el0g  t  log  k  - 


2«  n  , 
TToi¥  log  \t 


%t)]  eif?/ 


4)Sgn  t 


(A- 38) 
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— Continued, 


(A- 39) 


or 


u(t) 


"2  e  -  (log  t)  7  2  log  k 


logt 


(Sgn  t) 


1/2 


{A- 40) 


It  becomes  apparent  that  (A- 40)  is  only  an  approximation  when  one  attempts  to  verify 
the  functional  relationship 


d 

dt 


u  (t) 


C  u  (kt)  . 


(A- 41) 


The  left  hand  side  is 


u  *  (t) 


log  t  ±  ,  2*n  1 
log  k  J  log  k  j 


(A-  4  2) 
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A.  4  — Continued. 


and  the  rii,hi  hand  side  is 


u  (kt) 


u  (t) 

t 


(A-43) 


so  that 


Const  x  u  (kt> 


<A-44) 


unless  t  =  1. 
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